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SECTION  I 
INTRODUCTION 

Certain  physics  and  engineering  problems  require  the 
use  of  model  or  reference  atmospheres  to  represent  the 
environment.  Such  problems  include  determinations  of 
engine  and  aircraft  performance,  aerodynamic  characteristics, 
skin,  compartment  and  equipment  temperatures  under  transient 
climb  conditions  and  calculations  associated  with  the  vertical 
ascent  and  descent  of  missiles  or  munitions.  These  particular 
model  atmospheres  must  not  only  be  hydrodynamically  consistent 
in  the  vertical  and  homogeneous  in  composition  such  as  in  a 
standard  or  reference  atmosphere  but  also  must  be  applicable 
on  a  worldwide  basis.  A  more  complete  discussion  of  standard 
and  reference  atmospheres  and  some  definitions  of  terms  used 
in  this  report  are  given  in  Appendix  E,  The  most  common 
standard  atmosphere  is  the  1976  U.S.  Standard  Atmosphere 
(Reference  10).  However  it  only  represents  idealized, 
steady-state-conditions  near  45°  latitude,  primarily  over 
land  areas  and  does  not  accurately  represent  conditions 
at  any  given  place  or  at  any  given  time  of  year,  season, 
month  or  day. 

Similarly,  model  atmospheres  specifying  the  vertical 
envelopes  of  environmental  extremes  have  been  useful  in  the  design 
of  weapon  systems  to  operate  under  various  conditions  of 
temperature,  pressure,  density,  humidity,  wind,  etc.  The 
most  commonly  used  of  these  model  atmospheres  have  been 
the  "hot  day"  and  "cold  day"  atmospheres  taken  from  the  Hot 
and  Cold  Atmosphere  Tables  in  the  1957  MIL-STD-210A  (Refer¬ 
ence  6) .  These  tables  specified  the  temperature  extremes 
(hot  or  cold)  expected  at  various  altitudes  (levels)  from 
sea  level  to  100,000  feet.  However,  they  do  not  accurately 
estimate  conditions  likely  to  be  encountered  during  vertical 
motion  through  the  atmosphere.  These  atmospheres  were  actually 
constructed  to  provide  an  estimate  of  the  ten  percent  cal¬ 
culated  risk  for  the  hottest  and  for  the  coldest  areas  of 
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the  world,  level-by-level,  without  regard  to  the  re¬ 
lationship  between  levels.  Thus  the  vertical  temperature 
distribution  in  its  entirety  as  given  in  either  the  Hot  or 
Cold  Atmospheric  Tables  will  never  occur  at  any  given  time. 
For  example,  in  the  mid-latitudes,  both  hot  and  cold  ex¬ 
treme  temperatures  can  occur  at  the  same  time  at  different 
altitudes  over  the  same  location.  Thus  it  was  imperative 
(as  stated  in  MIL-STD-210A)  that  problems  dependent  on 
integrated  temperature,  pressure  and  density  over  an  alti¬ 
tude  range  be  solved  by  using  the  Polar  and  Tropical  Atmos¬ 
phere  Tables  which  were  also  given  in  MIL-STD-210A.  These 
latter  reference  atmospheres  were  both  homogeneous  and 
hydrodynamically  consistent  and  represented  average  January 
conditions  in  the  polar  regions  and  annual  mean  conditions 
in  the  tropics  respectively.  However,  these  did  not  satisfy 
the  design  engineer's  need  for  vertically  consistent  hot 
and  cold  extreme  temperature  profiles.  Thus  design  engineers 
continue  to  use  the  Hot  and  Cold  Atmospheres  from  MIL-STD- 
210A  as  well  as  associated  terminology,  i.e.,  "Hot  Day", 

"Cold  Day"  and  "Tropical  Day".  A  comparison  of  the  two  sets 
of  MIL-STD-210A  Atmospheric  Tables  as  depicted  in  Figure  1 
clearly  shows  how  unrealistic  the  Hot  and  Cold  Atmospheres 
are  when  taken  in  their  entirety. 

Subsequently  the  MIL-STD-210A  Hot  and  Cold  Atmosphere 
Tables  were  replaced  in  the  1973  MIL-STD-210B  (Reference  7) 
by  tables  of  the  1%,  5%,  10%  and  20%  high  and  low  temperature 
extremes  with  altitude  for  worldwide  operations.  These 
tables  likewise  do  not  represent  vertically  consistent 
profiles  of  the  atmosphere  and  are  strictly  envelopes  of 
level-by-level  extreme  conditions.  Thus  they  should  not 
be  used  to  consider  the  total  influence  of  the  atmosphere 
on  a  weapon  system  or  piece  of  equipment  during  its  trajec¬ 
tory.  For  such  applications  and  in  lieu  of  the  MIL-STD- 
210A  Polar  and  Tropical  Atmosphere  Tables,  MIL-STD-210B 
recommends  the  use  of  the  annual  tropical  atmosphere,  and 
winter  (January)  and  summer  (July)  atmospheres  for  each 
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15°  latitudinal  band  outside  the  tropics,  as  published 
in  the  "U.S.  Standard  Atmosphere  Supplements,  1966" 

(Reference  9)  .  Although  these  atmospheres  are  vertically 
consistent,  it  is  sometimes  difficult  to  select  which  one(s) 
to  use  in  a  particular  situation.  They  are  also  somewhat 
inconvenient  for  the  design  engineer  to  use  in  a  computer 
environment.  Moreover,  since  these  atmospheres  represent 
seasonal  mean  atmospheres  they  do  not  provide  truly  extreme 
hot  or  cold  day  temperature  profiles.  Such  latter  profiles 
normally  only  occur  for  a  few  hours  of  the  day  and  only 
affect  the  lowest  part  of  the  atmosphere  in  the  hottest 
(coldest)  areas  of  the  world.  Thus  neither  a  hot  day 
profile  in  the  tropics  nor  a  cold  day  profile  in  the  polar 
region  can  be  represented  by  means  of  data  taken  over  all 
hours  of  the  day  for  all  days  of  a  season.  However,  such 
"Hot  Day"  and  "Cold  Day"  profiles  are  of  use  to  weapon 
systems  design  engineers  as  is  evidenced  by  their  continued 
consideration  of  the  MIL-STD-210A  Hot  and  Cold  Atmosphere 
Tables. 

This  has  also  led  to  the  use  of  other  unrealistic 
model  atmospheres  such  as  the  constant-departure-from-stan- 
dard  vertical  temperature  profile  models  found  in  some  air¬ 
crew  flight  simulator  software.  Such  vertical  temperature 
models  are  not  only  improbable  but  also  become  unrealistic 
when  hot  or  cold  temperature  extremes  are  desired,  especially 
over  high  surface  terrain  heights.  Similar  problems  occur 
when  engineers  try  to  use  test  range  reference  atmospheres 
as  design  criteria.  Such  atmospheres  are  usually  constructed 
in  table  form  from  annual  or  monthly  means  of  data  taken 
at  standard  pressure  levels  throughout  all  hours  of  the  day. 

In  addition,  tables  of  extremes  at  various  levels  will  be 
provided.  For  certain  high  altitude  ranges,  table  values  may  also 
include  data  for  levels  below  the  height  of  the  range.  Such 
data  may  be  unrealistic,  particularly  if  it  was  generated 
bv  downward  linear  extrapolation  of  range  surface  data  values. 

Thus  range  reference  atmospheres  do  not  provide  the  engineer 
with  vertically  consistent  extreme  temperature  profiles. 
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As  a  result.,  a  computer  program  has  been  written  to  provide 
engineers  and  simulator  software  designers  with  a  simple  means 
to  generate  model  reference  atmospheres  in  which  the  above 
deficiencies  are  resolved.  Reference  atmospheres  can  now  be 
generated  for  various  geographic  latitudes  by  varying  the  input 
ground  level  temperature  and  sea  level  pressure  values.  These 
atmospheres  are  not  only  homogeneous  and  hydrodynamically 
consistent  but  also  can  incorporate  hot  or  cold  ground  temperature 
extremes  as  well  as  terrain  height.  The  algorithms,  assumptions 
and  meteorological  basis  in  the  development  of  this  program  are 
discussed  in  Section  II. 
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SECTION  II 
PROGRAM  DESCRIPTION 

Given  a  geometric  altitude,  a  terrain  height,  a  ground  level 
temperature  and  an  altimeter  setting  (sea  level  pressure  in  inches 
of  Hg) ,  the  AMAG  program  will  compute  a  corresponding  internally 
consistent  pressure  altitude  and  ambient  air  temperature.  From 
this  output  the  user  can  calculate  most  of  the  environmental 
parameters  needed  to  solve  particular  engineering  problems.  For 
computational  purposes  all  temperatures  are  considered  to  be 
virtual  temperatures.  An  extensive  list  of  such  parameters  and 
the  appropriate  equations  for  calculating  them  are  provided  in 
Appendix  B. 

This  program's  terrain  following  feature  requires  the  input 

of  a  terrain  height  and  a  ground  level  temperature.  However,  if 

a  terrain  height  of  zero  is  input,  the  ground  level  temperature 

then  becomes  a  sea  level  temperature.  When  the  standard  sea  level 

temperature  (15°C),  zero  terrain  height  and  standard  sea  level 

pressure  (29.92  In.  Hg)  are  the  inputs,  the  output  will  be  the  1976 

U.S.  Standard  Atmosphere  temperature  and  pressure  altitude  for 

any  given  geometric  altitude.  The  allowable  range  of  geometric 

altitude  is  from  -2.0  km  to  32.0  km.  The  program  considers  the 

difference  between  geometric  altitude  and  geopotential  altitude 

to  be  negligible  since  for  aeronautical  purposes  (altitudes  below 

33  km)  this  difference  is  very  small  (a  maximum  difference  of  0.5% 

at  33  km) .  This  program  assumes  that  the  air  behaves  as  a  perfect 

gas,  is  completely  homogeneously  mixed  and  is  in  hydrostatic 

equilibrium.  The  allowable  range  of  terrain  height  is  from  -2.0  km 

to  5.9  km.  The  allowable  range  of  input  ground  level  temperatures 

is  from  -50.0°C  to  60.0°C.  Mean  Sea  Level  (MSL)  corresponds  to  the 

effective  value  of  the  earth's  radius  (6356.766  km)  at  which  the 

2 

acceleration  of  qravity  equals  9.80665  m/s  and  where  the  geometric 
altitude  equals  zero  geopotential  km.  A  unique  feature  of  this 
program  is  that  both  terrain  and  geometric  heights  below  sea 
level  are  allowed  as  well  as  geometric  heights  below  the  terrain 
height.  This  feature  was  included  to  allow  independency  and 
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flexibility  in  selecting  input  geometric  and  terrain  height  values. 
For  example,  if  this  program  was  used  in  real-time  aircrew  flight 
simulator  software,  this  program  would  not  cause  the  main  program 
to  "error  off"  when  the  aircraft  geometric  height  was  below  the 
terrain  height. 

The  program  uses  a  terrain  following  five-layer  vertical 
temperature  model  based  on  the  1976  U.S.  Standard  Atmosphere 
and  the  1966  U.S.  Standard  Atmosphere  Supplements.  The  model 
contains  a  2  km  thick  ground  radiation  boundary  layer  and  a 
tropopause  whose  height  and  temperature  vary  according  to  the 
input  ground  level  temperature.  This  allows  the  user  to  specify 
very  hot  or  very  cold  ground  temperatures  typically  found  in  the 
tropics  or  arctic  regions  respectively.  The  vertical  temperature 
profile  is  constructed  from  the  ground  up  by  calculating  the 
temperatures  at  the  base  and  the  top  of  the  lowest  layer  and 
of  each  layer  in  sequence.  The  program  only  constructs  as  much 
of  the  vertical  profile  as  is  necessary  to  compute  the  temperature 
and  the  pressure  altitude  at  the  required  geometric  altitude. 

Within  each  layer,  a  constant  lapse  rate  of  temperature  is  assumed. 
Thus  the  temperature  at  a  given  height  can  be  calculated  either  by 
linear  interpolation  with  height  between  the  temperatures  for  the 
base  and  the  top  of  the  layer  or  by  upward  linear  extrapolation 
with  height  of  the  layer's  base  level  temperature  using  a  predeter¬ 
mined  vertical  temperature  lapse  rate. 

The  pressure  altitude  is  also  calculated  from  the  ground  up 
by  correcting  the  geometric  altitude  for  the  variation  of  the 
input  sea  level  pressure  (altimeter  setting)  and  for  the  deviation 
of  the  vertical  temperature  distribution  from  that  assumed  in  the 
standard  atmosphere.  For  the  latter  correction,  the  mean  tempera¬ 
ture  of  each  layer  is  compared  to  the  mean  temperature  of  that 
layer  in  the  standard  atmosphere.  The  corrected  thickness  values 
for  each  layer  are  then  added  together  to  arrive  at  the  pressure 
altitude. 
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SECTION  III 

PROGRAM  DESIGN  AND  SPECIFICATIONS 

The  construction  of  the  AMAG  model's  various  layers 
as  depicted  in  Figure  2  is  as  follows: 

a.  The  first  and  lowest  layer  is  the  terrain  layer 
which  is  bounded  at  the  bottom  by  -2.0  km  and  at  the  top 

by  the  input  terrain  height  (ground  level) .  The  temperature 
everywhere  within  this  layer  is  set  equal  to  the  input  ground 
level  temperature.  The  portion  of  the  tenain  layer  that 
is  below  sea  level  is  not  considered  in  the  pressure  alti¬ 
tude  computations  except  when  the  geometric  height  is  both 
below  sea  level  and  below  the  terrain  height. 

b.  The  second  layer  is  a  2  km  constant  thickness 
terrain  following  radiation  boundary  layer.  The  base  of 
this  layer  is  the  top  of  the  terrain  layer  at  which  the 
temperature  is  the  input  ground  level  temperature.  The 
method  of  calculating  the  temperature  at  the  top  of  this 
layer  is  rather  empirical  in  nature,  based  on  the  author's 
experience  and  on  the  need  to  account  for  both  hot  and  cold 
ground  level  temperature  extremes.  This  method  is  described 
as  follows: 

(1)  In  this  program  it  is  assumed  that  all  ground 
level  temperature  extremes  (those  less  than  0°C  or  greater 
than  30°C)  are  related  to  an  excess  or  deficit  of  surface 
heating  and  are  confined  to  the  lowest  2  km  above  ground 
level.  The  need  for  a  simple  linear  vertical  temperature 
lapse  rate  model  within  this  radiation  boundary  layer  led 
to  the  choice  of  2  km  as  its  thickness.  The  2  km  value, 
although  seemingly  high,  is  borne  out  in  the  model  reference 
atmospheres  contained  in  the  1966  U.S.  Standard  Atmosphere 
Supplements  (See  Figures  3,  4,  and  5).  Furthermore,  between 
2  km  and  8  km,  the  vertical  temperature  profiles  for  these 
latter  atmospheres  tend  to  depart  only  by  15°C  at  the  most 
from  the  1976  U.S.  Standard  Atmosphere  and  also  have 
comparable  tropospheric  vertical  temperature  lapse  rates. 
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ALTITUDE  (Geopotential  km  above  MSL) 


The  magnitude  and  algebraic  sign  of  this  departure  appear  to  be  a 
function  of  the  difference  between  the  sea  level  temperature  and 
the  standard  atmosphere  sea  level  temperature  of  15 °C.  Thus,  for 
this  model,  it  is  assumed  that  sea  level  temperatures  less  than 
0°C  or  greater  than  30°C  are  extremes  and  that  the  2  km  boundary 
layer  provides  a  simple  method  for  connecting  such  extremes  to 
the  2  km  height  where  the  vertical  temperature  lapse  rate  then 
becomes  equal  to  that  of  the  standard  atmosphere.  From  a  physical 
point  of  view,  2  km  is  the  upper  limit  of  surface  boundary  layer 
friction  and  diurnal  temperature  effects. 

(2)  To  relate  the  model's  ground  temperature  input  to  a 
sea  level  temperature  requires  the  calculation  of  an  "equivalent" 
sea  level  temperature.  When  the  terrain  height  is  zero  or  positive, 
this  temperature  equals  the  ground  level  temperature.  However,  if 
the  ground  level  temperature  is  greater  than  30 °C,  then  the 
"equivalent"  sea  level  temperature  is  set  equal  to  30°C.  Likewise, 
if  the  ground  level  temperature  is  less  than  0°C,  the  "equivalent" 
sea  level  temperature  is  set  equal  to  0°C.  When  the  terrain  height 
is  negative,  the  "equivalent"  sea  level  temperature  is  calculated 
by  extrapolating  the  ground  temperature  upward  to  sea  level  using 
the  standard  atmosphere  lapse  rate  and  then  correcting  as  necessary 
for  low  or  high  extremes  to  0°C  or  30 °C  respectively.  The 
temperature  at  the  top  of  the  2  km  boundary  layer  is  the  calculated 
by  extrapolating  the  "equivalent"  sea  level  temperature  upwards 
or  downwards  using  the  standard  atmosphere  lapse  rate.  Temperatures 
at  any  level  within  this  2  km  boundary  layer  can  be  calculated  by 
linear  interpolation  with  height  between  the  temperatures  at  the 
top  and  at  the  bottom  of  the  layer. 

c.  The  third  layer  is  bounded  at  the  bottom  by  the  top  of 
the  boundary  layer  and  at  the  top  by  the  model’s  tropopause. 

Within  this  tropospheric  layer,  the  temperature  lapse  rate  is 
always  equal  that  of  the  standard  atmosphere.  The  height  of  the 
top  of  this  layer  is  a  function  of  the  height  of  the  tropopause. 

Here  it  is  assumed  that  cold  sea  level  temperatures  are  associ¬ 
ated  with  a  low  altitude  warm  tropopause  and  warm  sea  level 
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temperatures  are  associated  with  a  high  altitude  cold 
tropopause.  Fmpirical  formulas  were  derived  to  calculate 
the  tropopause  height  and  temperature  directly  from  the 
model's  "equivalent"  sea  level  temperature.  These  formulas 
interpolate  linearly  between  those  values  for  the  tropopause 
heights  and  the  tropopause  temperatures  which  are  assumed 
to  be  associated  with  the  "equivalent"  sea  level  temperatures 
of  0°C,  15°C,  and  30°C.  Thus,  for  15°C  the  standard 
atmosphere  tropopause  height  and  temperature  values  of  11  km 
and  -56.5°C  are  assumed.  For  30 °C  the  values  used  are  16 
km  and  -74°C  and  for  0°C,  they  are  8  km  and  -52°C.  These 
specific  values  were  chosen  for  simplicity  and  to  mesh 
perfectly  with  their  coincident  "equivalent"  sea  level  temp¬ 
eratures  using  the  standard  atmosphere  lapse  rate.  These 
values  also  reasonably  fit  annual  averages  of  corresponding 
values  given  in  the  1966  U.S.  Standard  Atmosphere  Supplements 
and  facilitate  the  output  of  the  1976  U.S.  Standard  Atmosphere 
when  standard  sea  level  conditions  are  used  as  input.  The 
temperature  at  any  level  within  this  layer  can  be  calculated 
by  linear  interpolation  with  height  between  the  temperature 
values  computed  at  the  bottom  and  at  the  top  of  this  layer. 

d.  The  fourth  layer  is  bounded  at  the  bottom  by 
the  tropopause  and  at  the  top  by  a  fixed  height  of  20  km. 

This  layer  represents  the  lower  stratosphere  and  has  an  isother¬ 
mal  temperature  profile.  The  temperature  at  any  level  within 
this  layer  is  equal  to  the  tropopause  temperature  which  is 
calculated  from  the  "equivalent"  sea  level  temperature.  The 
latter  is  also  used  to  calculate  the  tropopause  height.  The 
value  of  20  km  as  the  height  of  the  top  of  this  layer  is  the 
same  as  that  used  in  the  1976  U.S.  Standard  Atmosphere. 

e.  The  fifth  layer  represents  the  upper  stratosphere 
and  is  12  km  thick,  extending  from  20  km  to  32  km.  The  layer's 
upper  boundary  height  and  temperature  values  are  the  same  as 
those  used  in  the  1976  U.S.  Standard  Atmosphere,  i.e.,  32  km  aM 
-44.5°C  respectively.  With  the  layer’s  lower  boundary  height 
fixed  at  20  km  and  its  temperature  equal  to  the  tropopause 
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temperature ,  the  temperature  at  any  level  within  this  layer  can 
be  calculated  by  linear  interpolation  with  height.  The 
temperature  lapse  rate  in  this  layer  is  thus  a  function  of  the 
tropopause  temperature  or  in  reality  a  function  of  the 
"equivalent"  sea  level  temperature. 

The  temperature  profiles  for  4  specific  cases  as 
calculated  by  the  AMAG  Program  are  depicted  in  Figures  6  and  7 . 
Specific  input  data  values  for  each  case  are  provided  in  the 
figure  legends.  It  should  be  noted  that  the  AMAG  program  does 
not  actually  calculate  complete  vertical  profiles  but  only 
computes  one  temperature  value  and  one  pressure  altitude  value 
corresponding  to  one  input  geometric  height.  A  vertical  profile 
of  temperature  and  pressure  altitude  could  be  generated  by 
placing  the  subroutine  AMAG  calling  statement  in  an  iterative 
loop  involving  geometric  height.  However,  one  should  be  very 
careful  about  interpolating  between  output  values  since  the 
input  geometric  heights  may  not  correspond  to  the  model 
atmosphere's  breakpoint  temperature  heights.  In  these  cases  it 
would  be  better  to  use  the  desired  geometric  heights  directly 
as  inputs. 

Some  caution  should  be  used  when  interpreting  the 
results  of  this  model  when  the  surface  temperature  inputs  are 
extremes  within  a  given  latitudinal  band.  This  model  (particularly 
the  tropopause  algorithm)  was  basically  designed  to  simulate 
the  mean  annual  temperature  profiles  (+  one  standard  deviation) 
for  various  latitude  bands  such  as  those  given  in  the  1966  U.S. 
Standard  Atmosphere  Supplements.  Thus,  for  example,  the  model 
cannot  always  simulate  representative  atmospheres  for  either  a 
very  cold  day  in  the  subtropics  or  for  a  very  hot  day  in  the 
arctic.  Likewise,  diurnal  temperature  effects  are  not  actually 
accounted  for.  However  the  model  will  produce  profiles  very 
similar  to  the  MIL-STD-210A  polar  and  tropical  atmospheres. 

See  Figure  8  for  a  comparison  of  this.  In  Figure  9  the  model 
profiles  are  also  compared  to  constant-departure-from-standard 
temperature  profiles. 
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ALTITUDE  (Geopotential  km  above  MSL) 


SECTION  IV 

PROGRAM  COMPUTER  CODE 

The  AMAG  program  has  been  written  as  a  subroutine  in 
Portrain  IV  language.  The  actual  program  listing  and  flow  chart 
are  given  in  Appendix  A.  The  program  uses  °C,  feet  and  In.  Hg 
as  units  of  temperature,  height,  and  pressure  respectively 
but  can  easily  be  modified  to  use  either  metric  or  English 
or  British  Engineering  units  exclusively. 

Execution  time  for  5000  calls  to  a  compiled  version  of 
this  subroutine  on  a  CDC  Cyber  175  is  approximately  one  second. 
Program  length  on  the  CDC  Cyber  175  is  313  60-bit  words. 

Questions  concerning  this  code  should  be  addressed  to  the 
ASD  Staff  Meteorology  Office  (ASD/WE) ,  Wright-Patterson  AFB,  OH 
45433. 


20 


SECTION  V 

PROGRAM  APPLICATIONS 


The  AMAG  program  output  parameters  of  pressure  altitude 
and  virtual  temperature  can  be  used  to  calculate  most  environ¬ 
mental  parameters  of  use  to  design  engineers.  The  equations 
necessary  to  do  this  are  given  in  Appendix  B. 

Additionally,  by  assuming  a  generalized  atmospheric 
vertical  moisture  model,  the  AMAG  output  virtual  temperatures 
can  be  converted  into  realistic  temperature  and  associated 
mixing  ratios  (absolute  humidity) .  The  appropriate  equations 
for  doing  this  are  given  in  Appendix  C. 

Furthermore,  Appendix  D  describes  a  synthetic  vertical 
atmosphere  scalar  wind  profile  model.  In  this  model  the  jet 
stream  height  is  made  to  be  consistent  with  the  AMAG  tropopause 
height  and  a  specified  surface  wind  input. 

This  AMAG  program  currently  is  limited  to  geometric  altitudes 
below  32  km.  With  minor  changes  to  the  program  computer  code,  it 
can  be  extended  to  47  km  geometric  altitude. 
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SECTION  VI 
CONCLUSIONS 

The  AMAG  program  has  been  tested  under  a  variety  of  input 
conditions.  The  results  show  that  AMAG  can  reproduce  exactly 
the  1976  U.S.  Standard  Atmosphere  and  good  fits  to  the  U.S. 
Standard  Atmosphere  Supplements,  1966  (MIL-STD-210B)  as  well 
as  the  MIL-STD-210A  Polar  and  Tropical  Atmospheres  and  several 
test  range  reference  atmospheres. 

Within  the  limitations  and  assumptions  required  to  develop 
a  simple  universal  program  of  practical  use  to  engineers,  the 
AMAG  program  can  be  successfully  applied  to  many  meteorological, 
physics  and  engineering  problems  including  simulation. 

Furthermore,  it  is  believed  that  this  program  can  bridge 
the  gap  existing  between  the  MIL-STD-210A  and  MIL-STD-210B 
model  atmospheres. 
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APPENDIX  A 
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PROGRAM  LISTING  AND  FLOW  CHART 
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SUBROUTINE  AMAG  (ZHGT*THGT»TTMP»ALTSTG»ZTMP,ZPA»iER> 

C  _ _ 

C******* 

C  _____  _  _ 

C  "  . . .  ”  ’  '  '  '  . . 

C _ PURPOSE -  - -  __ -  - - 

C  GIVEN  A  GEOMETRIC  HEIGHT  ♦  Z  H  G  T  ♦  I~N  FEET,"  A  “TERRAIN  HEIGHT  ♦fHGT* 

C  IN  FT, A  terrain  TEMPERATURE  *TTMP*  IN  OEGREES  C  AND  AN  ALTIMETER 

C  "  S FT T I NG  *ALTSTG*  IN  INCHES  CF  HG«  THIS  PROGRAM  USES  AN  EMP IR I CALLY 

C_  DERIVED  MODEL  ATMOSPHERE  TO  COMPUTE  THE  APPROXIMATE  TEMPERATURE 

C  . +ZT*<o*  in  OEGREES  C  ANC  THE  PRESSURE  ALTITUDE*  *ZPA*  IN  FEET.’  BOTH 

C  _  TERRAIN  A^NO  CECMETRIC  HEIGHTS  BELOW  SEA_LEVEL  AND  ALSO  GEOMETRIC 
C  '  -"heTghTS  PE-LOW  THE  T F R P A f Ni"  HE  I CHTS  “ARE  UNIQUELY  ALLOWED. 

C 

c  VARITblE  OFSCRfPMCN  --- 


c 

VARIABLE 

TYPE 

I/O/C 

DEFINITION  AND  PURPOSE  OF  PROGRAM  VARIABLE 

c 

IER 

INT. 

OUTPT 

EOUALS  0  IF  NO  INPUT  ERPGR,  CTHERWISF  ERROR 

c 

ALTSTG 

S.P. 

INPUT 

ALTIMETER  SFTTING  IN  INCHES  OF  HG 

c 

THGT 

S.P. 

INPUT 

TFPRAIN  HEIGHT  IN  FFET 

c 

ZHGT 

S.P. 

INPUT 

GEOMETRIC  HFIGHT  ABOVE  mfan  SEA  LEVEL  IN  FT 

c 

TTmo 

S.P. 

INPUT 

TERRAIN  TFHPEPATURE  IN  DEGREES  C 

c 

Z  P  A 

S.P. 

OUTPT 

PRESSURE  ALTITUDE  IN'  FEFT 

f 

ZTi-P 

S.P. 

OUTPT 

AMBIENT  TEMPERATURE  IN  OEGREES  C 

c 

W1 

S.P. 

CALC. 

STANCARO  ATMOSPHEPE  TEMPERATURF  IN  DFG  C 

c 

W2 

S.P. 

CALC. 

EQUIVALENT  SEA  LEVEL  TEMPFRATUPE  IN  DEG  C 

c 

W3 

S.P. 

CALC. 

GEOMETRIC  HEIGHT  ABOVE  E5616.fi  FEET  IN  FEET 

r 

W6 

S.P. 

CALC. 

MAXIMUM  HEIGHT  ABOVE  BOUNDARY  LAYER  IN  FEFT 

c 

W  5 

S.P. 

CALC. 

GEOMETRIC  HFIGHT  ABOVF  TERRAIN  HEIGHT  IN  FT 

c 

Wfc 

S.P. 

CALC. 

GEOMETRIC  HFIGHT  OF  TOP  OF  ROUNDARY  LAYER 

c 

W  7 

S.P. 

CALC. 

TROPOPAUSF  TRMPERATURF  IN  QFGREES  C 

c 

Vfi 

S.P  . 

CALC. 

ZTMP  in  DEGREES  K 

c 

W9 

S.P. 

CALC. 

standard  atmosphfre  tempepaturf  Wl  IN  DEG  K 

c 

W1C 

S.P. 

CALC. 

barometer  height  in  feet 

r 

Wll 

S.P. 

CALC. 

geometric  .fight  of  TRCPOPAUSF  in  feet 

c 

Wl? 

S.P. 

CALC. 

TERRAIN  TFmPFRATURE  TTMP  IN  DEGREES  K 

c 

wn 

S.P. 

CALC. 

STANDARD  ATMOSPHERE  TEMPERATURE  AT  TOP  OF 

c 

TERRAIN  LAYEP  IN  DEGREFS  K 

c 

VIA 

S.P. 

CALC. 

AVER  AGF  STANDARD  ATMOSPHFRF  TEMPERATURE  FOR 

c 

SEA  LEVEL  AND  TCP  DF  TERRAIN  LAYER 

c 

V  l  5 

S.P. 

CALC. 

FCUI VALENT  SEA  LEVEL  TEMPERATURE  W2  IN  DEGK 

c 

WIN 

S.P. 

CALC. 

STANCARO  ATMOSPHERE  TF>PF°ATURE'  AT  TOP  OF 

r 

BOUNDARY  LAYER  IN  CEGRFCS  K 

C 

W  1  7 

S.P. 

CALC. 

TFMPFRATURE  AT  TOP  OF  BCUNOARY  LAYFP.  IN 

r 

CFGRFFS  w 

c 

wie 

S.P. 

CALC. 

TROPOPAUSF  TFMPERATURF  W7  IN  DEGREES  K 

r 

*19 

S.P. 

CALC. 

MAXIMUM  CF  GEOMETRIC  HCT  VS  TERRAIN  HGT 

C 

W?C 

S.P. 

CALC. 

MINIMUM  OF  W5  VS  2  KM 

C 

V?1 

S.P. 

CALC. 

MINIMUM  OF  ZHGT  VS  Wll 

C 

WZ2 

S.P. 

CALC. 

“IN!  mu'"  OF'  ZHGT  VS  STD  AT>OS  TROPOPAUSF  HGT 

C 

C  INPUT  AND  TUTPUT  REQUIREMENTS  - 

C  INPUT  THPCUCH  ARGLMFNT  LIST 

C  "  OUTPUT  THROUGH  ARCUmENT  LIST 

c 

r  rfStr  icTYcns'--- . 

C  ZHGT  "UST  LIE  BETWEEN  -6560.  AND  106987.  FEFT 

C  T H G T ’ H 0 S T  LIF  BFTWEFN  -6560.  AND  19357.  FEET 
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C  TT“P  •list  LtF  BETWEEN- -SO.'O  AND  60.0  DEGREES  C*  . 

f  ALTSTC  “(JST  LIE  BETWEEN  28.00  AND  31.00  INCHES  OF  HG 

C 

C  ERROR  CCNOTTIHAS  AND  RETURNS  - 

C  ”  IE>-1  IF  7HGT  IS  0UTSI0F  ALLOWABLE'RANGE 
C  IFR«?  IF  THOT  IS  OUTS  I  OF  ALLOWABLE  RANGE 

f  iER -A  IF  T  TCP  IS  OUTS  I OE  ALLOWABLE  RANGE-  "  . . 

C  IFR*fi  IF  ALTSTC  IS  OUTSIOE  ALLOWABLE  RANGE 

C  IE  IFR>F,  PROCEAM  RETURNS  IMMEDIATELY  WITH  IEP»SUM  CF  IER  VALUES" 

C 

C  ACCURACY  - 

C  TF-PCPATUPC  IS  ACCURATE  TO  WITHIN  1  OFGPFE  C  A^O  PR ESSURE_ALT_I_TUOE 

C  Ttj  W  i TH I m”  ICO  FEET  . ”  ~ . . . . 

C 

C  RFFFRFNCFS - 

C  U.S.  STANDARD  AT*<CSPHFRF,  1076 

f  19*6  U.S.  STANTARO  ATHOSPWFRE  SUPPLEMENTS 

C  A  S  0- TP- 79- S  c  S (  — AUTOMATED  MOCEL  AT^OSPHFRF  generator  PROGRAM 

c  '  (  A  **  A  G  )  —'ref  1979  '  . .  . . 

r 

C  LANGUAGE .  . . 

C  r f P T P A N  IV 

C  . . 

C  CnRF  STGPAGF  -  LESS  THAN  320  WORDS  DECIMAL 

c  * .  '  . .  "  . . .  . . . 

r  CIS*.  CR'J"  OR  TAFT  PPOUIRFBENTS - NOME 

C 

C  EXECUTION  T I  **F - 1.0  SFC  FOR  5C00  CALLS 

C 

r  PPFRFnuiSITr  PPOCPA«S/SUPR0UTINES - A*AX1  AND  AMIN1  FUNCTIONS 

C  "ACHINF  rFPFAr/FNCE - NON'F 

c  "  ' 

C  PRECISION  -  SINCLF 

r 

r  AOCITIfNAl  FNTPY  POINTS  -  NONF 

r 

r  ***************************************************************  ******** 

r 

C  CUTPUT  r  A  T  A  INITIALIZATION 
0 


C 

C 

r 


c 

r 

r 

c 

C 


7 p A  »  O.C 
2T«P  *  C.C 
IFR  *  0 


INPUT  PATA  f^ROP  CHECK 

IF  1 7*-O-T.LT.-*E60..CR.7EGT.GT. 104937.  >  IFP  *  IFR  ♦  I 
IE  »THGT.LT.‘-f*A.n..aR.TEGT",GT.iqV>7.»  IER  »  IER'  ♦  2 
IF  »TT*'P.LT.-EC..OP  .TTwP.GT.fcC.  )  IEP  *  IER  ♦  B 
IF  (ALTSTC.LT.2B.PC.QR.ALTSTr.GT.Tl.00)  IFR  >  IEP  ♦  0 
!c  (IFR)  10.K,99 


CALCULATE  STANDARD  AT-OSP^ERF  TEMPERATURE  W1  USING  the  1tD__ATM0SPHERE 
LAPSE  p  A  T r  "E  C.CCL0B12  OECPFFS  C  AFOOT,  T«0P0PA'LSF"  TFMPFfcAfURE  OF 
-■if.'  OTF  0  ♦  TRTPPPAUSF  HEIGHT  CE  36089,237  FEET 


mis  pass  is  msi  »u»Mn  i SAsnwai 

Q(>Pf  FUS8USKB5D  TO  90C 


t 


oono  r»  r»  ^  non  o  r>  o  o  o  o  r*>  o  o  r>  o  c~>  o 


r,  CAICULATE  HEIGHT  ABOVE  65616.  796  FEET 

C  _  _ 

10  W3  *  A*A*1 (C.C  ,7HCT-65616. 7R61  . ~ 

Hl_  *  15.  -  AHIM(7i.5.0.0010812*ZHGT)  ♦  ( C._Q0030A8*W3 1 _ 

C  "  . 

C  CALCULATE  PAPQ»ETER  HE IGHT_T9  HEIGHT  OF  STANDARD  ATMOSPHERE  SEA  LEVEL 

c  pressure  surface'  ' '  . . . . . .  . . . 

C  _  _ _ 

H10  *  030.“*  (20.02  -  ALTSTG 1 

check ”if~s7Tnoarc  atmosphere  temperature" oesTrecT. 

fp“T  rp'o  f  ;  F  6 . 0  .’  C  .  A  So  .  T T  M  p  .  E  C  .  15.  7  ~G  fi  “T  0  IV 

CALCULATE  EQUIVALENT  SEA  LEVEL  TEMPERATURE  W2  (SE E"«ASlc'REF?PFNCE )  " 

w?  ft-p  •  .  . . .  . 

IJF  (THGT.LT.C  ,CI  W2  «  W?  ♦  <  0 . 00 1  OBI  2  ♦THGT  ) 

If  TtmcT:lt.c'.c  .ano.  w2.FQ.i'V.o  '.and.  ThgtIgeTthgtI  or' 'to  15 

IF  (W2.LT.0.0I  W2  »  C.O 

IF  YvP.rXlC'.C  »'  W2  »  3C.0  * . 

GO  TC  20  _  . . 

FOP  STANDARD  AT^CSPHFPE.  AGO  JURCWETRIC  HEIGHT  to  GEOMETRIC  HEIGHT  TO 

get  Pressure"  "altitude;  then'V  return. . . 

15  If  HP  *  W1 

7P A  «  7HGT  ♦  MO 
Rc  TURN 

CALCULATF  CrOM'FT'R'lC  HEIGHT  ABOVE  GROUND  LEVEL  W5  . . 

2C  W5  »  7HCT"  -  THCT 

CALCULATF'  GEOMETRIC  HFIGHT  OF  TCP  OF  BOUNDARY  LAYEP 

~  WE '*'  THC'T"*'f5tr.fcB  . .  .  . 

CALCULATF  TFMPEPATURf  FOR  CEQ“ETR  IC  HEIGHT  BELOW  TERRMNXfVEL" 


IF  (V5.CT.0. »  GO  TO  26 
2 T"P  •  TTMp 

~  WTO  "3? . .  ”  '  *  " . ”  . . . . . 

CALCULATE  MAXIMUM  OF  INPUT  HEIGHT  VS  TOP  OF  6561.68  FOOT  BOUNDARY  LYR 
26.  WA  «  A  •*  A  Y  l  (  7HC  T  »  W6  ) 

fALCULATF~YTwPfRATUPE  FOR  IN^UT  "HEIGHT"  ABOVE  BOOnCARY  IaYER'JXE'HP 
AT  THF  TROPPPALSE  ANO  HEIGHT  OF  THF  TRDPOPAUSE 

IF( (W2-16.I.GT.0.)  GO  TC  20 

W7  '*  -5  2.0  ( C’.  3  ♦  W  2  1 

Wll  *  262*6,718  ♦  ( 656.  16706*W2 1 

"7TMP  -  AK*xi  (fc7,W2-r0.0Cr08l2*W*)T  ♦  W  3>T  if.  00010 0"VVc7oF<Fo7  6  2  ♦  W2  ) 

GO  TC  3C 

20  W7  «  -30‘,C  -  n>2'*7./6.)  ~  . . .  . . 
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r 

c 

c 

c 

r 

C 

r 

c 

r 

r 

C 

r 

c 

c 

c 

r 

C 

r 

C 

r 

r 

C 

r 
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C 

C 

C 
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c 
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wu  *  lSEPfi.02e  ♦  (  lO<n.Hl32*W2) 

7T“P  *  A»Ayi»V7,W2-IC.0C19812*WA) )  ♦  W3*(-0.0001397»0.00002963*W2) 

CALCULATE  Tf*PFPATU»F  WHFN  WITHIN  BOUNDARY  LAYER  BY  INTERPOLATION 
BETWEEN  VAl.UcS  AT  THE  0  AS  E  AND  TOP  QE  TUE  BOUNDARY  LAYER 

30  Ic  ( /^OT.LT.WE  1  71*9  *  TTMp  ♦  U W5 /6S6 1 . 68  I  * < 7 TMP-T T^pTI ~ 

CCNVEP  T  TEMPER  A TLCCS  TO  DECREES  K 

Wifi  **»?♦  > 73, IE  '  “ 

37  we  «  i'»P  ♦  273. IE 
W9  a'Wl  ♦  273.1  5 
Wl?  »  Tfwp  ♦  77'».3E 
US  *  w?  *  ?7?. IE 

CALCULATE  STtNCARC  ATMOSPHERE  TEMPERATURE  AT  TOP  OF  TERRAIN  LAYER  ' 
wii •;*>»> Vis  -  (c.ooi9ei2*TWGT) 

CALCULATE  AVERA.r.f  STANDARD  ATw0SPHF»F  TEMPERATURE  BETWEEN  SEA  LVL  AND 
TCP  OF  TERRAIN  LAYFP 

VI*  *  (2Pfi.lS  ♦  W 1 3  >  /  2  _  . . 

CALCULATr  STANCACr  ATNCSPUFRF  TEMPERATURE  AT  TOP  OF  POUNDARY  LAYER 
HA  »  M3  -  13. 

CALCOLATF  Tr*pf.o.M(jOF  AT  TOP  DF  ROUND  ARY  LAYER 
V 1 7  *  VIS  -  <r.C019fil2AW6) 

CALCULATE  PPTSSURF  ALTITLOF  7PA  RY  ADDING  PARO**ETER  HE  I CHT  To  THE 
SUW  OF  TFT  PRESSLPE  ALTITUDE  THICKNESSFS  FOR  EACH  SUCCESSIVE  VERTICAL 
LAYFP 


7PA 


W10 


CALCULATE  THE  PPESSURF  ALT  I TUOF  THICKNESS  FOR  A  GIVEN  LAYER  BY  MULTI¬ 
PLYING  the  LAYFR  CFO“EtqiC  THICKNESS  BY  THE  RATIC  WITHIN  THE  LAYER  OF 
THE  AVFRACE  STANTAPO  A  T  MQS  P  HE0  F  TEMPERATURE  TO  THE  AVERAGE  AMBIENT^ 

tevpfratuRf 

CHECK  IE  TFPRAIN  HFIGHT  ANC  GEQNFTPIC  HEIGHT  BOTH  ABGVE  SFA" LEVEL 
AC  IF  (THCT .CT  .G.C.AND.7HGT .GT .0.0)  GO  TO  60 

CHECK  IF  TFPRAIN  HEKHT  AT  OP  BFLEIW  SEA  LEVEL  ANt  IF  'g60>TTRU  HEIGHT 
AROVF  SEA  LEVEL 

a  1  IF  UHCT.Lf  .C.C.ANC.7HGT.GT.O.O)  GO  TC  55 

CHECK  ^jr  TFPRAIN  HEIGHT  AND  GFOmeTPIC  HEIGHT  BOTH  AT  OR  BELOW  SEA  LVL 
A2  IF  (THGT.Lf  .C.C. AND. 7HGT.LE. 0.01  CO  TC  A8 
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ym*  cop*  furkisked  k>  wxj 


C  SINCE  "TFPRAJN  HEIGHT  IS  A8CVF  SFA  LFVFL  AND  geometric 'height  is  at  or 
C  PELCW  SFA  LEVEL,  CALCULATE  ZPA  AND  RETURN 
f 

A3  ZPA  »  7P4  ♦  <  ZHGTK  28 R. 1 K- ( C. 0009906*7HCT )) /W 12 ) 

'return . 

C  _  _  _  _  _  _  _  _ 

c  'CALCULATE  ZPA  FOR  POUNGAPY  LAYER  BELOW  SEA  LEVF’L 
C 

A?  W19  «  A P A y l ( 7FCT,THOT) 

Z 0 A  «  ZPA  ♦  W19*(?88.1S-0.0C099C6*W19)  / 

£  ~  "lVil2~+WP—THGT*(W17-W  121/656 1.68)  ♦  2 

C  _ 

C  CHECK  IF  CFP*'F TR 1 C  HEIGHT  IS  PFLCW  THE  TERRAIN  HEIGHT  WITH  BOTH  THE 
C  TERRAIN  HEIGHT  ANC  GEOMETRIC  HEIGHT  AT  OR  RFLCW  SEA  LEVFL 
C 

IF  (V5.CT.0.)  GC  TO  51 
C 

C  CALCULATE  ZPA  FOR  TERRAIN  L AYF  R  ANO  RETURN 

c  . . . 

ZPA  •  ZPA  ♦  n>«5*(wn-(O.OCO«J<»06*W5>  )/W12' 

51 'RETURN 

c 

C  SINCE'THF  TFRPAIN  HEIGHT  IS  AT  OP  BELOW  SEA  LEVFL  ANO  THE  GEOMETRIC 
C  HEIGHT  IS  AP.TVR  f>CTH  SEA  LFVFL  ANO  THE  TERRAIN,  PRE-CORRECT  FOR  TERR 
C  HEIGHT 'P FlCW _SFA  LEVEL  ANC  JUMP  TO  THE  BOUNDARY  L A YF R  CALCULATION 
C 

/PA  .  ZPA  ♦  THGT 
GO  TO  ZC 
C 

C  CALCULATE  ZPA  FHR  TERRAIN  LAYER  WHEN  both  GEOMETRIC  HFICHT  ANO  TERR 
C  "'HEIGHT  *AR*>  AB'CVF'SFA  level' 

c 

Vo'/P'A  *  >PA*APIM  C  THGT,  ZHGT  )*(  288.  15-0.000R906AAK  INK  THGT,  ZHGT)  )/Wl2 

c  _ 

c'  CHFCK  I F  GfT'-FTRIC  HEIGHT  IS  A80VF  HFIGHT  OF  TERRAIN 
C _  _ 

'  "if  iw*.lf;'o.  for  turn 

c  _ 

C  CALCULATE  ZPA  FOR  BOUNDARY  L A Y F P  WHF N  GEOMETRIC  HFIGHT  Is  ABOVE" SEA 
C  LH  ♦  A a  V  T  HF  TERRAIN  HFIGHT 
C  ‘ 

7C  W20  •  AMN1IWF, 6*61. 681 

*  ZP  A  .'“tr"A  ♦  W<C*(W13-O.C0099Q6YW?0)/(WlZ+w20*(Wl7-Wl2)/13i?3.36) 

c 

C  ‘CHECK  IF  GFC»F  TR If  HFIGHT  ABOVE  BOUNDARY  LAYER 
C 

7*  IF  KZHCT-WM.LF.O.  )  RETURN 
C 

C  C  At  CUL  AT  F  ZPA  FOR  TROPOSPHERIC  LAYER  FROM  THE  TO?  OF  THE  BOUNDARY’  LYR 
C  TO  THF  LCVEP  OF  FITHFR  The  STANOAPO  ATMOSPHERE  TROPOpAUSF  CR  THE 

C  TROPCPALTF 

C 

C  CHECK  IF  TPOPOPALSF  HEIGHT  IS  LESS  THAN  STANDARD  ATMOSPHERE  TROPOPAUS 
C 

"  IE"  livif-T6CPS.237)  ,  GT  .0.1  GC  TO  TO  . . 

c 

C  CALCULATE  ZPA  FOR  TRCPCSPHFRE  BEV.FEN  THf-  BOUNDARY  LAYER  ANO  TROPOPAUS 
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‘  C .  ■  . "  '  ' 

wn  >  Jl'IMI  ZHCT »W1 1 ) 

‘  ZPA  PA  i  (V21-W6)*(Wl~6N288.i5-C.  0019812  *W21 )  . 

£  /IW174-WIS-0.0C19812*W21» 

C 

C  _CHECK_IF  GEOMETRIC  WEIGHT  IS  RFLCW  TROPOPAUSE 

C  . .  . . 

81  IF  H7HCT-Wll).LE.O.)  RETURN 
C 

C  CALCULATE  ZPA  F®  C  M  TPPPOPAL'SF  TC  STANDARD  ATMOSPHERE  tropopause_ 

C  .  ■ 

ZPA  «  ZPA  M ( AM  IN  1 <  ZHGT  »38089.237)-Wll)*(288. 15-1 C. C 009906+ 1 W lit 

£  A  M~f  N 1  f /HOT  »3AiC89.237)  )  )  )  /  W 1  8  )  ~~ .  * . 

C 

c"  CHECK"  I F  GPOKETR'IC  HEICHT  is  ABOVE  STANDARD  ATMOSPHERE  iROPOPAUSE 

c 

8 A  IF  <<7Hr,T-3F089.237>.LE.0.>  RETURN 

0 

f"  “CALCULATE"  7~PA‘~FGP  LAVER  aETWEEN  ”$ T  AND AR D  ATMOsBhEpT  TRCpCPAUSE'  AND  " 
C  85tl£.79A-  FFET 

c‘  '  . . . . 

ZPA  •  ZPA  ♦  (  AMJN1I  7HCT-3f  089. 237, 29927. 558)*21A>j>A>5/W18)  _ 

8  7  GC  ™C  9f. 

r 

C  ''CALCULA"-»E"7PA  FCR  tROPOSPHFp'F  'fiFTvTFF'N  THE  "TOP  OF  THE  "BOUNDARY  LAYER  ♦ 
C  THE  STANDARD  AT“CSPHFPE  TRCPDPAUSF 

r 

90  W?2  -  AM»M  (  ZHfT,  16089.777) 

•'  ZPA"  i  'ZPA  U»22-W£)*<W16*?88.l5-C.O0  198 1 2*W2i?T . .  '  " 

£_  _  /<W17*W15-0.0019812*W2?> 

C  CWFOK  IF  OFOMF  TR  l  C  HEIGHT  IS  RELOW  STANDARD  ATMOSPHFRE  TROPOPAUSE 

91  IF  (I  ZHCT-3A-OF9.237J.LE.O.  )  RFTURN 

c  ‘  '  . .  '* 

C  CALCULATE  z PA  cop  LAYER  RETWFEN  STANCARD  ATMOSPHERE  TROPOPAUSE  AND 

C  TRtWPAUSF . .  .  . . . 

C 

ZPA'.  ZPA  ♦(  (A»TM1(  ZHGT,  Will  -  3’6089  . 2 1 7 1  *2 1 <E  .f  5  /  (Vi  5—”  " 

£  (C.CCr.99C6*ntCP9.237»AMlNi  (7HGT,Wll)  >1) » 

C  . 

f  CFFCK_JC  Oc^xciPIC  HEIGHT  IS  AROVF  TROPOPAUSE 

IF  ( I 7UGT-W 1 1  ).LF .0.1  R  F  TUP  N 

C  . . 

C  CALCULATE  /PA  fop  layer  °E7WFFN  TROPOPAUSE  ANO  N5616.796  FEET 

c .  . . 

7  R  A  •  ZPA  Ml  AMNII  ZHGT,  6S6 1  6 . 79f> ) -W  11  )  ♦?  1  <• .  68_/W  18) 

C 

C  CHFCK  IF  GE  P VF  TP  J  C  HFIGHT  IS  AROVE  98918.796  FFET 

C.  . 

96  IF  frfl.tr.O.)  RETURN 

C  '  '  . 

r  CALCULATE  ZPA  pr.o  LAYFP  ARCVE  65616.796  feet  _ _ 

ZPA  »  ZPA  ♦  (V3*(?1A-.E8*W9)/IW18NW8I) 

99  re  tip* 

END 
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Calculate  standard 
atmosphere 
temperature 

1 1  r 

Calculate  hgt  above 
20km  (lf<0,  set=0) 


Calculate 
tropopause 
height  and 
temperature 

Calculate 
temperature 
for  higher  of 
either  boundary 
layer  top  or 
(geometric  height 


(!) 


Calculate 

.  YES  ^ 

pressure 

- j 

> — — w- - 

altitude 

Return  ^ 


Set 

temperature 
equal  to 
terrain 
temperature 


Figure  A-l. 


Program  Flow  Chart 
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/  geometric''  v  N0 
.height  above  topo>— ) 
>vof  boundary/' 
'sj.ayer?/' 


Calculate 
temperature 
for  within 
boundary  layer 


Convert 
temperatures 
to  decrees  K 

Start 

pressure  altituds 
calculat ions 


Calculate 
standard 
atmosphere 
temperature  for 
boundary  layer. 

Set  totaf  pressure 
altitude  equal  to 
barometer  height. 

^Ar^v 

terraiiS^ 
weight  and  geometric 

N  height  above  / 
^ea  level>^ 


Add  terrain  layer 
pressure  altitude 
thickness  to 
total  pressure 
altitude. _ 


/terrainV 
/height  or  aK 
'below  sea  level 
&  is  geometric 
v  height  above  . 
Nsea  level>'' 


Pre-correct 
total  pressure 
altitude  for  portion 
of  boundary  layer 
below  sea  level. 


/'terrain  N. 

'and  geometric\N° 
heights  both 
at  or  below  / 
Nsea  level’/ 


_ z - - - 

Add  pressure  altitude 

thickness  for  below 
sea  level  boundary 
layer  to  total 
pressure  altitude. 


Since  terrain 
height  is  above 
sea  level  and 
geometric  height 
is  below  sea 
level,  add  pres¬ 
sure  altitude 
thickness  for 
terrain  layer  tc 
total  pressure 
altitude 

C  Return  ^ 


■''geometries, 
height  aboveX.  NO 
the  terrain .z'"” 
\he  i  gh  t 


''geometries, 
height  below' 
the  terrain > 
S^heightTx^ 


Return 


Add  boundary 
layer  pressure 
altitude  thickness  to 
total  pressure  altitude. 


Add  terrain  layer 
pressure  altitude 
thickness  to  total 
pressure  altitude. 


Return 


Figure  A-l.  Continued 
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c 


^  geometric\N 
height  above 
boundary  y 
""v.  layery/^ 

Tves 


Return 


./tropopauseN. 
^•''height  less  thanX^ 
^  that  of  standard 

atmosphere  y/ 
^Xfr  opopaus  e£/ 

YES 

Xdd  pressure  altitude 
thickness  for  layer 
between  boundary  layer 
and  tropopause  to 
total  pressure 
altitude 


Add  pressure  altitude 
thickness  for  layer 
between  boundary  layer 
and  standard  atmosphere 
tropopause  to  total 
pressure  altitude. 


/ geometries, 
height  above" 
tropopause^ 
Njieight/'^ 


Return 


>/geometriX. 

^  height  aboveS. 
standard  atmospherj 
S.  tropopause? 


Return 


.  ■  ■  . 1- 

Add  pressure  altitude 

Add  pressure  altitude 

thickness  for  layer 

thickness  for  layer 

between  tropopause  and 

between  standard  atmo- 

standard  atmosphere 

sphere  tropopause  and 

tropopause 

- c— - 

tropopause  to  total 

pressure  altitude. 

geometries. 

^  height  aboveS. 
standard  atmosphere 
's^tropopause^/ 


Return 


/geometries^ 
height  above 
ytropepause/ 


Return 


Add  pressure  altitude 
thickness  for  layer 
between  standard 
atmosphere  tropopause 
and  20km  to  total 
pressure  altitude. 


Add  pressure  altitude 
thickness  for  layer 
between  tropopause  and 
20km  to  total  pressure 

_ altitude 


Fiaure  A-l. 


Continued 


|  Add  pressure  altitudq 
thickness  for  layer 
between  20km  and 
geometric  height 
to  total 

pressure  altitude 


Return 


Figure  A-l.  Concluded 


APPENDIX  B 

PROGRAM  OUTPUT  APPLICATIONS 


Using  the  AMAG  subroutine  input  values  of  altimeter  setting  ALSTG 
(In. Hg),  geopotential  height  h  (feet),  terrain  height  h  (feet) 

and  ground  virtual  temperature  T^  (°C)  and  output  values  of 

pressure  altitude  hp  (feet)  and  virtual  temperature  T^  (°C), 

several  environmental  parameters  *  of  use  to  design  engineers 
can  be  calculated  as  follows: 


Standard  Atmosphere  Sea  Level  Value  of 
Temperature  T  (°K) 


=  288.15 


Pressure  P 


(lb/sq.  foot) 


=  2116.217 


Density  p  (lb/cu.  foot) 


=  0.076474 


Speed  of  Sound  a 


(feet/second)  =  1116.45 


Coefficients  of 

Dynamic  Viscosity JJ0  (lb/foot* second)  =  1.2024213  X  10  5 
Kinematic  Viscosity^  (sq.  feet/second)  =  0.1572284  X  10  ^ 
Thermal  Conductivity  k,  (BTU/foot • second • ®R)  =  0.4067468 

L  /  O 

Ambient  Temperature  T^  (°K)**  =  T^  (°C)  +  273.15 

Standard  Atmosphere  Temperature  T^  (°K)**  = 

Pst 

288 . 15-MIN (71. 5, 0.00 19812 *h  ) +0.000304 8*MAX (0.0, h  -65616 . 796) 

P  P 

*See  Table  B-l  for  conversion  factors  for  metric  to  English  units. 
**In  the  AMAG  model,  TK  —  T^  ,  PK  =  PK  and ^  by  definition. 


T,  0  T,  unless  h=h  .  However,  T.  normally  is  not  a 

p  .  nst  P  nst 

1  st 


meaningful  parameter. 
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w  * 

I 


t, . 


'  ^  **  X  > 


**  Jiv/* 


Ambient  Pressure  (lb/sq.  foot)** 
For  hp  ±.  36089.237 


=  P0  fT, 
Qst 


or  PQ  / 

stt 


145442-hr 


-5.25585 


5.25585 


145442 

For  36)389.237  <  hp  <  65616.796 
=  472 . 678exp  (1.73456-0.0000480631  h_) 


For  hp  > 


65616.796 
=  114.343  f  216.65/Th 

P 


0 


st 


or  114.343 


t 


710794. 


645177.2  +  hr 


^  34. 

-) 


163 


34.163 


** 


=  0.010413414  jph/Th(°K)J 
=  Th  (°K)/288.15 


=  Ph/211C.217 


.  76474 


Ambient  Density  (lbs/cu.  foot) 

Ambient  Temperature  Ratio  T^/Tq 

f\  s  t 

Ambient  Pressure  Ratio 

Ambient  Density  Ratio 

Pressure  Altitude  Variation  PAV 
(feet) 

Speed  of  Sound  a^  (feet/second) 

Mach  Number  M  for  a  given  speed  V  in  =  v/e?^ 

feet/second  or  0.0152044  v/  £t(°K)|  1//2 

Density  Altitude  h^  (feet) 

For  h  £  36089.237 


■  hp-h 


=  65.77 


[VKj| 


1/2 


=  145442 


[l-<%**Z0*3*9] 


For  36089. 237  <  h  <  65616.796 
P 


=  36089.237-20807.0  log. 


For  h  >  65616.796 
P 


V0. 0227188  / 


=  710794  /0.0054981V'0284389  -645177.2 

{  Sh  1 


For  h  <  10000  feet, 
P 


h  + 
P 


120  [Th  (°K)  -  Th 


(°K)] 


st 
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il 


Total  Temperature  Ratio  Tfc/X 
Free  Stream  Total  Pressure  to 


=  1.0  +  0.2M" 


Ambient  Pressure  Ratio  (Subsonic) Pc  =  £l . 0+0 . 2M^| 3 ’ 


Total  Pressure/Ambient  Pressure 
Ratio  Pt  (Supersonic)  from 
~P 

Rayleigh  Pitot  Equation 
Incompressible  Dynamic  Pressure 

(lb/sq.  foot) 


=  166.92M2/[?M2-1.0]  5/2 


Compressible  Dynamic  Pressure  t* 

(lb/sq.  foot) 

v 

Equivalent  Air  Speed  V  (feet/second)  =  V 


=  0.7PhM 

or  1481.352 
or  0.0001618 3 /phV 

tar 

or  0.5  V2 

=  p|^(l.0+0.2M2)3*5-1.0] 


or 


Calibrated  Air  Speed  V  (feet/second) 


29.00751^ 


/2 


For  Subsonic  Flow 


=  2496.646 


For  Supersonic  Flow 
=  984.66 


-1.0 


1/2 


C+  v) 


5/2 


1/2 


u  '  -*  /  \ 

Coefficient  of  Thermal  Conductivity  ^  ^ *  \  h/ 

(RTIJ/f oot •  second  •  0 T.  +  245.4  X 


3/2 


(BTU/foot* second • °R) 

Coefficient  of  Dynamic  Viscosity 
( lb/foot • second) 


Coefficient  of  Kinematic  Viscosity  72. 
(sq.  feet/second)  * 


xirrTO7^T~ 

4  4 

/j.97973X10~6(Th)  3/2 
Th+110.4 
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TABLE  B-l 

CONVERSION  FACTORS 

To  Convert  From  To 

Length 

Meter  foot 

U.S.  statute  mile 
U.S.  nautical  mile 


Mass 

Gram 


pound  -  mass 
slug 


Volume 

Cubic  meter  cubic  inch 
cubic  foot 


Pressure 

Millibar  Pascal  (Newton/sq.  meter) 

millimeter  Hg. 
inches  Hg. 
pound/sq.  inch 
pound/sg.  foot 

Speed 

Meter/second  foot/second 

kilometer /hour 

mile/hour 

knot 


Density 

Gram/cubic 

meter 

Temperature 

"Celsius  (C) 


Coefficient  of 

Kinematic  Vis- 
cosity(Sq. 

meter/sec^  sq.  foot/sec 

Coefficient  of 

Dynamic  Vis¬ 
cosity  (Newton • 

second/sq.m)  Ib/foot’ second 

Coefficient  of  Thermal 
Conductivity 

(Watts/meter  •  9K)  BTU/foof  second  •  °R 


pound/cu.  foot 
slug/cu.  foot 


"Fahrenheit  (F) 
"Kelvin  (K) 
"Rankine  (R) 


* 


Multiply  By 

3.2808398  +0* 
0.0621371  -2 
0.0539957  -2 


2.2046224  -3 
0.0685217  -3 


6.1023759  +4 
3.5314667  +1 


1.0000000  +2 
7.5006151  -1 
2.9529971  -2 
1.4503768  -2 
2.0885437  +0 


3.2808398  +0 
3.5999997  +0 
2.2369363  +0 
1.9438447  +0 


0.0624279  -3 
0.0019403  -3 


°F  =  1.8C  +32 
"K  =  C+273.15 
"R  =  1.8K 


1.0763909  +1 


0.6719689  +0 


0.1606044  +3 


NOTE: 


Indicates  the  power  of  10  for  positioning  the  decimal 
point . 


APPENDIX  C 

ATMOSPHERIC  MOISTURE  MODEL 

In  this  Appendix,  a  method  is  outlined  for  generating  vertical 
atmosphere  moisture  profiles  which  are  consistent  with  any  AMAG 
output  virtual  temperature  -  pressure  altitude  profiles. 

The  method  involves  a  three-layer  vertical  moisture  model 
based  on  an  average  relative  humidity  vertical  profile  derived 
from  the  U.S.  Standard  Atmosphere  Supplements,  1966  (see  Table  E-l) . 
It  uses  the  AMAG  output  values  of  virtual  temperature  and  pressure 
altitude  and  the  AMAG  input  values  of  ground  temperature,  geometric 
height,  terrain  height  and  altimeter  setting. 

The  first  and  lowest  layer  is  the  tropospheric  layer  bounded 
by  mean  sea  level  and  the  tropopause.  For  heights  at  and  below 
mean  sea  level  a  constant  80%  relative  humidity  is  assumed.  Above 
mean  sea  level  the  relative  humidity  is  assumed  to  exponentially 
decrease  with  height.  The  terrain  height  is  purposely  ignored  as  a 
level  of  moisture  discontinuity  since  surface  moisture  is  not 
usually  conserved  over  large  changes  in  terrain  height.  High 
elevation  locations  tend  to  be  much  drier  than  low  elevation 
locations.  In  this  layer  relative  humidity  is  used  as  the  moisture 
modeling  parameter,  rather  than  a  parameter  which  specifies  a  certain 
amount  of  moisture  such  as  specific  humidity.  This  is  done  to  avoid 
the  possibility  of  encountering  supersaturated  conditions  in 
calculations  involving  cold  tropospheric  temperatures.  Relative 
humidity  (RH)  at  any  positive  geometric  altitude  h  for  below 
the  tropopause  is  calculated  as  follows: 

RH  (%}  =  80.  exp  (-h/8  km)  (C-l) 

The  tropopause  height  h^  is  calculated  by  the  AMAG  algorithm 
involving  ground  temperature  as  follows: 

For  ground  temperatures  T^  less  than  0°C, 

h  =  26246.718  feet  =  8  km  (02) 

For  ground  temperatures  between  0°C  and  15°C  (T  in  °C)t 

h.  =  26246.718  +  (656.16796  *  T  )  feet  (C-3) 
t  g 

For  ground  temperatures  between  15°C  and  30?C  (T  ir.  0  C), 

h  =  19685.038  +  (1093.6132  *  T  )  feet  (04) 
x.  y 

For  ground  temperatues  T^  above  30°C^ 

ht  =  52493.434  feet  =  16  km  (C-5) 
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Since  the  tropopause  tends  to  act  as  a  cap  on  the  upward 
propagation  of  moisture  from  the  troposphere,  the  atmosphere 
above  the  tropopause,  i.e.,  the  stratosphere, is  very  dry.  The 
discontinuity  in  moisture  between  the  troposphere  and  stratosphere 
occurs  at  what  is  called  the  hygropause  which  is  located  at  a 
height  of  one  km  above  the  tropopause.  Above  the  hygropause  and 
up  to  100,000  feet,  the  amount  of  moisture  is  extremely  small, 
i.e.,  a  mixing  ratio  of  one  to  three  parts  water  vapor  mass  per 
million  parts  of  dry  air.  Thus  in  this  model  we  will  assume  a 
constant  mixing  ratio  at  all  altitudes  above  the  hygropause. 

From  this  assumption,  the  second  and  third  layers  of  the  moisture 
model  can  be  defined. 

The  second  moisture  layer  lies  between  the  tropopause  and  the 

height  of  one  km  above  the  tropopause.  The  moisture  modeling 

parameter  used  in  this  layer  is  the  mixing  ratio  which  is  defined 

as  the  dimensionless  ratio  of  the  mass  of  water  vapor  to  the  mass 

of  dry  air.  The  term,  saturation  mixing  ratio,  is  that  value  of 

the  mixing  ratio  which  corresponds  to  completely  saturated  air, 

i.e.,  100%  relative  humidity.  Thus  relative  humidity  is  defined 

as  the  percent  ratio  of  the  mixing  ratio  to  the  saturation  mixing 

ratio.  The  mixing  ratio  at  the  tropopause  is  calculated  from  the 

tropopause  temperature  Tfc,  tropopause  height  hfc,  tropopause  relative 

humidity  RHfc  and  tropopause  pressure  altitude  PAfc.  The  tropopause 

temperature  T.  is  calculated  from  the  ground  temperature  T  as 
u  y 

follows : 


For  ground  temperatures  T^  less  than  0°C, 

Tfc  =  -52  °C  =  221 . 15°K 

For  ground  temperatures  between  0°C  and  15°C  (T  in  °C) , 

y 


(06) 


T,  =  -52  °C  -  (T  *  0.3)  (07) 

t  g 

For  ground  temperatures  between  15°C  and  30°C  (T  in  °C) , 

T  =  -39 °C  -  (T  *  7./6J  (03) 

t  9 

For  ground  temperatures  above  30°C/ 

T  =  -74  °C  =  199 . 15°K  (09) 

The  tropopause  height  hfc  and  relative  humidity  RHt  are  calculated 
using  the  formulas  given  previously  in  this  appendix.  The  tropo¬ 
pause  pressure  altitude  PA^.  can  be  calculated  by  the  AMAG  program 
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using  the  tropopause  height  h^.  as  the  input  geometric  altitude  along 
with  the  ground  temperature,  terrain  height  and  altimeter  setting 
inputs.  From  PAt,  the  tropopause  pressure  Pfc  can  be  calculated 
using  the  formulas  given  in  Appendix  B.  At  this  point  the 
tropopause  mixing  ratio  rfc  can  now  be  calculated  by 


t  f)* 


{±11 £ _ \ 

V  Tt  +  237  y 


RHt  X  10 


(C-10> 


where  P  =  tropopause  pressure  in  mb 


RHt  =  tropopause  relative  humidity  in  percent 
T  =  tropopause  temperature  in  °C 
The  mixing  ratio  at  the  hygropause  and  for  the  moisture  model '  s 
third  layer  can  now  be  determined .  It  is  set  equal  to  the 
smaller  of  either  the  tropopause  mixing  ratio  rfc  or  the  value  of 
0.003  gm/kg.  For  the  first  case,  the  mixing  ratio  at  all  geometric 
altitudes  above  the  tropopause  is  set  equal  to  the  tropopause 
mixing  ratio.  For  the  second  case,  the  mixing  ratio  r^  for 
geometric  altitudes  between  the  tropopause  and  the  hygropause  is 
calculated  by  logarithmic  interpolation  with  height  between  the 
value  of  rfc  and  of  0.003  gm/kg  as  follows: 

In  r^  (gm/kg)  =  (ht(km)-  h(km)+  lkm]*£ln  rt(gm/k£)+  5 . 8 0 8 8 J -5 . 8088 

(C-ll) 

Several  relationships  exist  between  relative  humidity,  mixing 
ratio,  virtual  temperature  and  temperature.  For  geometric  height 
below  the  tropopause,  the  temperature  T  can  be  estimated  from  the 
AMAG  output  virtual  temperature  Tv  by 


(°c)  a*  t 


(°C)  -  10* 


fTv-14)/38] 


(C-12) 


This  approximation  is  necessary  for  T  values  above  -25°C.  For 
values  of  Tv  less  than  -25°C,  the  second  term  is  negligible 
and  TfcTv*  The  approximation  avoids  solving  a  complex  exponential 
equation  and  allows  the  direct  calculation  of  mixing  ratio  rh  for 
any  given  height  h  from  the  temperature,  pressure  altitude  and 
relative  humidity  by  the  following  equation: 
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a*  ' 


■■^aawattaaaBaflesflaesftiigaa^ 


(£g)#38*0 


RHh  *  10 


(.  ’-1  ’»  ) 

\T.  +  237.3/ 


A 


(C-13)* 


where  RH^  =  relative  humidity  in  percent 
=  temperature  in  °C 

=  pressure  in  mb  as  calculated  from  the 
pressure  altitude  using  equations  given 
in  Appendix  B 

Once  r.  is  known,  then  a  more  accurate  value  of  T,  can  be 
h  h 

calculated  from  T  . 


V°K)=T 


<°K) 


[l  - 


(C-14) 


In  addition,  the  mixing  ratio  as  approximated  above  is  equal  to 
the  specific  humidity  which  is  defined  as  the  ratio  of  the  mass 
of  water  vapor  to  the  mass  of  moist  air  containing  the  vapor. 
Absolute  humidity  which  is  simply  the  density  of  the  water  vapor 
can  be  calculated  as  follows: 


A  fc)  - flD  *A  <8  * 10 

where  r^  =  mixing  ratio 


-3 


(C-15) 


A  =  dry  air  density  as  computed  from  equations 
for  ambient  density  given  in  Appendix  B. 
Example  moisture  profiles  are  provided  in  Table  C-l. 

The  tropospheric  moisture  profile  calculated  by  this  model  can  be 
made  drier  or  wetter  simply  by  adjusting  the  mean  sea  level  80% 
relative  humidity  value  to  lower  or  higher  values  respectively. 
However,  it  is  not  recommended  that  values  of  less  than  30%  be  used. 


*Equations  C-10  and  C-13  were  derived  from  Tetens'  empirical 
formula  for  saturation  vapor  pressure  over  water.  In  this 
appendix  all  formulas  and  calculations  were  based  on  saturation 
over  water.  To  make  calculations  valid  for  saturation  over  ice, 
replace  the  constants,  7.5  and  237.3,  by  9.5  and  265.5  respectively. 
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TABLE  C-l 
MOISTURE  PROFILES 


Geopotential 
Altitude  Above 
M.S.L.  (km) 

COLD 

DAY* 

HOT 

DAY* 

Absolute 

Humidity 

(gm/m^) 

Virtual 

Temperature 

(°C) 

Absolute 

Humidity 

(grn/rn^) 

Virtual 

Temperature 

(°C) 

0 

2.66-1-0** 

-5.0 

22.47+0** 

33.0 

1 

1.76+0 

-9.0 

15.39+0 

25.0 

2 

1.15+0 

-13.0 

9.47+0 

17.0 

3 

6.05-1 

-19.5 

6.09+0 

10.5 

4 

3.08-1 

-26.0 

3.81+0 

4.0 

5 

1.51-1 

-32.5 

2.21+0 

-2.5 

6 

7  ' 

7.16-2 

3.21-2 

-39.0 

-45.5 

1.42+0 

7.83-1 

-9.0 

-15.5 

8 

1.38-2 

-52.0 

4.34-1 

-22.0 

9 

1.34-3 

-52.0 

2.32-1 

-28.5 

10 

1.14-3 

-52.0 

1.19-1 

-35.0 

11 

9.76-4 

-52.0 

5.89-2 

-41.5 

12 

8.65-4 

-52.0 

2.78-2 

-48.0 

13 

7.17-4 

-52.0 

1.24-2 

-54.5 

14 

6.14-4 

-52.0 

5.26-3 

-61.0 

15 

5.26-4 

-52.0 

2.10-3 

-67.5 

16 

4.51-4 

-52.0 

7.81-4 

-74.0 

17 

3.86-4 

-52.0 

5.00-4 

-74.0 

20 

2.43-4 

-52.0 

2.99-4 

-74.0 

25 

1.11-4 

-48.88 

1.22-4 

-62.12 

30 

5.16-5 

-45.75 

5.29-5 

-50.25 

32 

3.81-5 

-44.50 

3.82-5 

-44.50 

Zero  terrain  height  and  29.92  altimeter  setting  is  assumed. 
Indicates  the  power  of  10  for  positioning  the  decimal  point 


APPENDIX  D 

ATMOSPHERIC  WIND  MODEL 

In  this  appendix,  a  method  is  outlined  for  generating 
vertical  atmospheric  wind  profiles  which  are  consistent  with 
any  AMAG  output* virtual  temperature-pressure  altitude  profiles. 

The  method  involves  a  six-layer  model  of  an  idealized 
steady-state  scalar  wind  speed  profile.  It  uses  the  AMAG  input 
values  of  geometric  altitude,  terrain  height  and  ground  temperature 
and  a  user-specified  value  of  the  ground  wind  speed.  This  wind 
model  is  terrain  following  such  that  below  the  terrain  height, 
the  wind  speed  is  set  to  zero  and  that  above  the  terrain  height, 
there  is  a  one-km  thick  boundary  layer  in  which  the  wind  speed 
increases  exponentially  with  height  to  a  value  of  one  and  one- 
half  times  the  ground  wind  speed  at  the  top  of  the  boundary  layer. 
Above  the  boundary  layer,  the  wind  speed  increases  with  height  to 
a  maximum  at  the  height  of  the  jet  stream.  The  height  of  the 
maximum  wind  or  jet  stream  is  a  linear  function  of  the  height 
of  the  tropopause  and  varies  from  8  km  to  14  km.  The  maximum 
wind  speed  is  set  equal  to  six  times  the  ground  wind  speed.  Above 
the  jet  stream  the  wind  speed  decreases  with  height  up  to  20  km 
where  a  wind  speed  value  of  twice  the  ground  wind  speed  is  reached. 
Between  20  km  and  23  km,  the  wind  speed  is  held  constant.  Above 
23  km  the  wind  speed  increases  linearly  at  a  fixed  rate  of  4m/sec 
per  kilometer. 

Wind  speed  can  thus  be  calculated  for  any  geometric  height 
as  a  function  of  the  ground  wind  speed,  terrain  height  and 
tropopause  height.  The  tropopause  height  is  calculated  from  the 
ground  temperature  as  in  Appendix  C. 

Let  =  wind  speed  at  geometric  altitude  h 

V  =  ground  wind  speed 
g 

Vj  =  maximum  wind  (jet)  speed  =  6V^ 

h^  =  tropopause  height 

h  =  terrain  (ground)  height 
9 

hT  =  jet  stream  height 
J 
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(D-l) 


then  for  h  <  h 


vh  =  0 


for  hg  £  h  ^(hg  +  lkm) 

Vh  =  Vg  e>:p  [  ‘ 4  *  (h_hg)  km”3J 
for  +  1  km)  <  h  <  hj 


h  =  0 .  75h  +  2  km 

J  1 


(D-2) 


(D-3) 


V.  =  6V 
h  g 


for  h  <  h  <20  km 

U 


V.  =  6V 
h  g 


for  20  km<  h  <23  km 

V.  =  2V 

h  g 

for  23  km  <  h  <•  32  km 


(hj  -  V  i) 


Lfr.j-ho-1)2  +  3  (Vh) 


r  < 

20 

km  -  hj  2 

L 

[20km- 

|2  +  2(h-hj)2_ 

V,  =  2V  +4 
h  g 


*(h'23  km) 


(D-4) 


(D-5) 


( D~  6 ) 


(D-7) 


Example  wind  profiles  are  provided  in  Figure  D-l. 

Since  this  wind  model  represents  only  idealized  steady- 
state  scalar  wind  conditions,  typical  of  mid  latitude  locations, 
its  application  should  be  limited  to  preliminary  design  investiga¬ 
tions.  The  wind  profiles  generated  with  this  model  are  approximately 
90  percentile  envelopes  without  considering  gust  factors,  wind 
shear  or  wind  direction.  For  more  advanced  design  work  for  specific 
operational  capabilities  at  specific  locations,  vector  synthetic 
wind  profiles  based  on  detailed  data  from  such  locations  should  be 
used  (See  Ref.  3,  Chap.  8).  If  a  wind  direction  is  needed  for 
convenience  sake,  it  should  be  assumed  to  be  towards  the  east  as 
is  typical  of  the  mid  latitude  westerlies,  i.e.,  a  wind  direction 
of  210° .  This  value  of  270°  can  vary  slightly  from  240°  to  300° 
if  northerly  or  southerly  wind  components  are  required.  The  wind 
direction  variation  with  height  cannot  normally  be  generalized  but 
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in  the  mid  latitude  westerlies,  wind  direction  tends  to  turn 
clockwise  with  height  between  the  boundary  layer  and  the  height 
of  the  maximum  wind.  Lastly  if  a  wind  profile  with  a  given 
maximum  wind  (jet)  value  is  desired  then  the  input  ground  wind 
speed  should  be  set  equal  to  one-sixth  of  that  value. 
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APPENDIX  E 

STANDARD  ATMOSPHERES* 


1 .  INTRODUCTION 

Numerous  problems  in  physics  and  engineering  are  sensitive  to 
the  state  of  the  atmosphere.  For  example,  aerodynamic  problems 
are  density  and  temperature  dependent;  refraction  problems  are 
dependent  on  density  and  water  vapor  concentration;  and  radiation 
transport  problems  are  very  sensitive  to  the  concentrations  of 
molecular  species  in  the  atmosphere. 

In  order  to  allow  investigators  to  normalize  the  solutions  of 
these  and  other  problems  to  common  atmospheric  conditions,  various 
standard  atmospheres  have  been  developed.  In  general,  they  represent 
an  idealized  model  of  the  mean  of  a  large  number  of  atmospheric 
measurements.  As  our  ability  to  accurately  measure  the  atmospheric 
parameters  of  interest  has  developed,  the  standards  have  changed 
somewhat.  Since  a  large  number  of  measurements  of  the  lower  few 
kilometers  of  the  atmosphere  have  been  available  for  many  years, 
standards  there  have  changed  little,  if  any,  for  the  last  20  years 
or  so.  However,  as  radiosondes  have  been  improved  and  other 
techniques  have  become  available,  measurements  have  been  made  to 
higher  and  higher  altitudes.  It  is  at  these  higher  altitudes  that 
most  of  the  recent  changes  have  been  seen. 

2.  DEFINITIONS 

At  this  point,  it  is  necessary  that  we  develop  a  common 
vocabulary  for  the  discussions  to  follow.  The  following  are 
definitions  commonly  accepted  by  atmospheric  scientists: 

a.  Standard  Atmosphere:  A  hypothetical  vertical  distribution 
of  atmospheric  temperature,  pressure,  and  density  which,  by 
national  or  international  agreement,  is  taken  to  be  representative 
of  the  atmosphere  for  mean  annual  conditions  at  45°M  latitude. 

The  air  is  assumed  to  obey  the  perfect  gas  law  and  the  hydrostatic 
equation.  It  is  further  assumed  that  the  air  is  dry  and  that  the 
acceleration  of  gravi'y  does  not  change  with  height. 

*This  Appendix  is  an  updated  extra ction  from  prog ram  documentation 
for  the  FTD  Standard  Atmospheric  written  by  Capt  John  D.  Mill  of 
the  FTD  Staff  Meteorology  Office  (FTD/WE)  in  1975. 
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b.  Model  Atmosphere:  Any  theoretical  representation  of  the 
atmosphere.  A  model  atmosphere  does  not  have  the  force  of  law  or 
international  agreement. 

c.  Reference  Atmosphere:  A  model  atmosphere  generally  accepted 
as  representing  the  atmosphere  under  certain  conditions  such  as  a 
particular  season,  or  latitude  or  altitude  interval.  They  are  often 
moist  atmospheres  and  may  include  such  things  as  concentrations  of 
molecular  or  atomic  species.  Many  are  intended  to  supplement 
atmospheres  and  therefore  often  have  been  called  supplementary 
atmospheres,  a  term  whose  definition  has  now  changed  to  refer  to 
such  things  as  ozone  or  humidity  distributions. 

d.  Temperature  (Kinetic  Temperature) :  A  measure  of  the  mean 
kinetic  energy  of  the  translational  motion  of  molecules  or  atoms. 

e.  Molecular-scale  Temperature;  A  fictional  temperature 
derived  from  application  of  the  perfect  gas  law  under  the  assumption 
that  the  mean  molecular  weight  of  dry  air  is  a  constant.  Below 
about  90  km,  it  is  equal  to  the  kinetic  temperature. 

f.  Virtual  Temperature:  A  fictional  temperature  of  moist 
air  derived  from  the  application  of  the  perfect  gas  law  under  the 
assumption  that  the  mean  molecular  weight  is  that  of  dry  air. 

This  temperature  is  often  used  in  reference  atmospheres  to  facili¬ 
tate  calculations  of  density  and  density-dependent  factors. 

g.  Mean  Molecular  "Weight”:  The  weighted  average  of  the 
molecular  masses  of  the  atmospheric  constituents,  excluding  water 

vapor  (dry  air) .  It  has  been  defined  on  the  carbon  12  scale  of 

.  4  -3 

atomic  mass  as  2.89644  x  10  kg*m  . 

h.  Geopotential  Height:  The  height  of  a  given  point  in  the 
atmosphere,  relative  to  sea  level,  proportional  to  the  potential 
energy  of  a  unit  mass  at  that  height.  It  arises  from  the  assump¬ 
tion  of  constant  acceleration  of  gravity  (see  paragraph  1) .  It 
is  defined  as: 
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where:  H  =  geopotential  height 

gg  =  acceleration  of  gravity  at  sea  level 

z  =  geometric  height 

g  =  acceleration  of  gravity  at  z. 

(The  standard  for  gn  varies  somewhat,  depending  on  the  model,  but 
.  u  -2 

is  most  often  taken  as  9.80665  m»sec  at  45°N  latitude.) 

i.  Altimeter  Setting:  The  value  of  atmosphere  pressure  to 
which  the  scale  of  a  pressure  altimeter  is  set.  This  setting 
represents  the  pressure  required  to  make  the  altimeter  indicate 
zero  altitude  at  mean  sea  level. 

j.  Pressure  Altitude:  The  altitude,  in  the  standard 
atmosphere,  at  which  a  given  pressure  will  be  observed.  It  is  the 
indicated  altitude  of  a  pressure  altimeter  at  an  altimeter  setting 
of  29.92  inches  of  mercury.  A  pressure  altimeter  converts  atmos¬ 
phere  pressure  into  altitude  using  standard  atmosphere  pressure- 
height  relations. 

3 .  RESULTS 

a.  Proper  Application  of  Standard  Atmospheres: 

Standard  atmospheres  are  idealized  representatives  of  mean 
conditions  near  45aN  latitude,  primarily  over  land  areas.  As  such, 
they  do  not  accurately  represent  conditions  at  a  given  place  or  at 
a  given  time.  The  primary  value  of  a  standard  atmosphere  lies  in 
its  use  as  a  reference  point  by  which  different  calculations  can 
be  compared.  The  adoption  of  a  standard  assures  that  differences 
in  results  are  due  to  elements  of  the  experiment  other  than  the 
atmospheric  parameters.  For  this  type  of  application,  it  does 
not  matter  a  great  deal  what  the  exact  standard  is,  but  only  that 
the  same  standard  is  used  for  all  calculations  which  are  to  be 
compared. 

On  the  other  hand,  there  are  many  problems  where  atmospheric 
conditions  are  of  central  importance.  Experimental  data  often  need 
to  be  reduced  to  common  environmental  conditions  (often  a  standard 
atmosphere) ,  or  the  sensitivity  of  some  value  to  changing  conditions 
may  need  to  be  determined  (often  relative  to  a  standard  atmosphere) . 
In  these  cases  —  and  others  —  actual  data  at  the  time  of  the 
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experiment,  specially  developed  model  atmospheres,  or  one  or 
more  reference  atmospheres  are  needed. 

The  degree  of  sophistication  required  in  environmental 
data  depends  on  the  experiment  (or  numerical  simulation)  in  a 
rather  complicated  manner.  The  basic  point  to  be  made  is  that 
a  careful  sensitivity  analysis,  error  analysis,  or  theoretical 
investigation  must  be  made  in  each  case.  It  may,  indeed, 
develop  that  a  standard  atmosphere  is  adequate,  but  this  must  be 
determined  —  not  assumed.  It  is  often  helpful  to  consult  an 
environmental  specialist  when  attempting  to  answer  this  question. 

The  point  of  contact  within  ASD  is  the  Staff  Meteorology  Office 
(WE) .  It  is  recommended  that  any  proposed  use  of  standard 
atmospheres  or  other  environmental  data  be  discussed  with  the 
staff  meteorologist. 

b.  Available  Standard  and  Reference  Atmospheres: 

In  this  section,  the  more  widely  used  standard  and  reference 
atmospheres  will  be  discussed  briefly.  There  are  a  number  of 
organizations  which  publish  standard  and  reference  atmospheres,  and 
there  is  a  rather  complex  interrelationship  among  them.  Many  of 
the  same  people  sit  on  committees  of  more  than  one  organization. 

The  abbreviations  used  for  atmospheres  of  these  organiza¬ 
tions  and  an  attempt  to  unscramble  these  relationships  follow: 

ARDC  -  Air  Research  and  Development  Command  (now  Air 
Force  Systems  Command).  The  earlier  U.S. 

Standard  Atmospheres  were  developed  and 
published  by  ARDC. 

ICAO  -  International  Civil  Aviation  Organization.  The 
ICAO  does  not  develop  its  own  atmospheres,  but 
adopts  those  developed  by  other  organizations, 
primarily  COESA. 

COESA  -  U.S.  Committee  on  Extension  to  the  Standard 

Atmosphere.  COESA  develops  recommendations  for 
revisions  to  U.S.  standards  and  has  been  very 
influential  in  their  adoption  as  international 
standards . 
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ISO 


-  In' ernational  Standards  Organization.  Coordinates 
the  work  of  various  national  organizations,  such 
as  COESA,  and  adopts  international  standards. 

CIRA  -  COSPAR* International  Reference  Atmospheres. 

COSPAR  does  not  develop,  recommend,  or  adopt 
standard  atmospheres,  but  publishes  a  number  of 
reference  atmospheres,  which,  in  general,  begin 
at  25  to  30  kilometers  and  extend  upward. 

CIRA  1965  -  Reference  atmospheres,  including  latitudinal 

(10°  intervals)  and  monthly  variations  from  25 
to  80  kilometers  and  mean  profiles  for  25  to 
500  kilometers  and  110  to  2000  kilometers. 

U.S.  Standard  Atmosphere  Supplements,  1966  -  (Hereafter 
referred  to  as  the  1966  Supplements.)  Moist  reference  atmospheres 
for  winter  and  summer  and  five  latitudes  (15°N,  30°N,  45°N,  60°N, 
and  75°N) ,  extending  to  120  kilometers  except  at  70°N  (30  km) . 
Includes  models  of  warm  and  cold  stratospheric  winter  regimes  for 
60°N  (to  80  km)  and  75°N  (to  30  km) ,  and  models  from  120  to  1000 
kilometers  for  various  levels  of  solar-geomagnetic  activity.  It 
is  very  similar  to  the  CIRA  1965,  except  at  high  latitudes,  where 
the  CiRA  1965  appears  to  be  in  better  agreement  with  recent  data. 

CIRA  1972  -  Revision  of  CIRA  1965.  The  latitudinal  and 
monthly  atmospheres  were  extended  to  120  kilometers. 

ISO  1972  -  The  latest  international  standard.  It  is 
identical  to  the  1962  Standard  (U.S.)  to  50  kilometers. 

U.S.  Standard  Atmosphere,  1976  -  Recently  published,  it 
is  identical  to  the  ISO  1972  to  80  kilometers  and  revises  the  1962 
Standard  above  50  kilometers. 

JACCHIA  -  A  series  of  models  extending  above  120  kilometers, 
developed  by  Dr.  L.  Jacchia  of  the  Smithsonian  Astrophysical 
Observatory.  Rather  than  tables  of  discrete  reference  atmospheres, 
they  are  algorithms  for  calculating  models  for  various  solar- 
geophysical  conditions.  In  their  various  forms,  they  provided 
the  basis  of  thermospheric  reference  atmospheres  since  1965. 

(CIRA  1965  and  1972,  1966  Supplements.) 

international  Council  of  Scientific  Unions  Committee  on  Space 
Research 
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Figure  E-l  represents,  graphically,  the  relationships  among 
the  major  standard  atmospheres  since  1952.  The  reference  atmos¬ 
pheres  are  too  numerous  to  include  here,  but  some  comparisons  may 
be  found  in  the  U.S.  Standard  Atmosphere  Supplements,  1966. 

The  following  is  a  representative  sample  of  the  large 
number  of  models  which  have  been  developed  in  the  past.  There 
were  standards  prior  to  1952,  but  they  were  essentially  identical 
with  more  recent  standards  at  low  levels  and  are  no  longer  used 
in  their  original  form.  In  all  cases,  the  more  current  edition 
replaced  the  older,  which  is  no  longer  considered  valid,  though 
often  it  is  unchanged.* 

ICAO  1952  -  A  standard  of  temperature,  pressure,  and 
density  to  11  kilometers.  It  has  been  extended  by  later  models, 
but  is  unchanged*  in  that  region. 

ICAO  1954  -  Extension  of  the  ICAO  1952  to  20  kilometers. 

It  is  also  unchanged*  bv  subsequent  standards. 

ARDC  1956  -  A  U.S,  standard  which  extended  the  ICAO  to 
500  kilometers. 

COESA  1958  -  A  proposed  revision  of  the  ARDC  1956  to 
300  kilometers.  It  was  never  adopted  as  a  standard  and  was 
superceded  by  COESA  1962. 

ARDC  1959  -  A  revision  of  the  ARDC  1956  above  53 

kilometers . 

CIRA  1961  -  The  first  COSPAR  reference  atmosphere.  It 
is  not  a  standard,  and  differs  somewhat  from  the  ICAO,  ARDC,  and 
COESA  standards. 

COESA  1962  -  Proposed  revision  of  the  ARDC  1959  above 
20  kilometers.  It  was  adopted  as  the  U.S.  Standard  Atmosphere, 

1962. 

U.S.  Standard  Atmosphere,  1962  -  (Hereafter  referred 
to  as  the  1962  Standard.)  A  revision  of  the  ARDC  1959  above  20 
kilometers,  it  includes  a  proposed  extension  of  the  ICAO  to  32 
kilometers.  It  also  includes  models  from  120  to  700  kilometers 
for  various  levels  of  solar-geophysical  activity. 

*  A 1  though  th  ere  haveHBeen  slight  revisions  in  some  physical  constants , 
the  basic  data  is  unchanged.  The  major  difference  is  due  to  the 
redefinition  of  0  C  as  273.15°K  vice  273.16°K.  In  this  paragraph, 
"unchanged"  is  given  this  meaning. 
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c.  Sources  of  Error 

As  pointed  out  a  number  of  times  above,  most  standard 
atmospheres  are  identical  in  the  lower  atmosphere.  In  general, 
all  standards  are  identical  up  to  20  kilometers.  When  choosing 
a  standard  for  use  above  these  altitudes,  it  is  important  to 
recognize  that  a  number  of  differences  exist  in  this  region, 
and  care  must  be  exercised  if  the  results  of  calculations  are  to 
be  compared  with  those  done  at  different  times  or  by  different 
investigators.  Where  standards  have  been  labeled  as  identical 
in  this  paper,  one  may  assume  that,  conceptually,  they  will  give 
identical  results  to  at  least  three  significant  figures.  There 
are  minor  differences  such  as  the  redefinition  of  the  ice  point  on 
the  Kelvin  scale  or  the  slightly  different  heights  of  some  reference 
points,  such  as  the  stratopause. 

Other  sources  of  error  involve  the  different  definitions 
of  height  and  temperature  given  in  Paragraph  2  above.  In  many 
tables  of  standard  atmospheres,  heights  are  given  in  both  geometric 
and  geopotential  meters  (or  feet)  and  care  must  be  exercised  when 
comparing  results.  The  difference  between  these  heights  varies 
with  both  latitude  and  altitude  and  varies  from  zero  at  sea  level 
to  several  kilometers  in  the  thermosphere  (above  120  km) .  For 
aeronautical  purposes,  the  maximum  error  can  be  taken  as  110 
meters  at  20  kilometers  near  the  equator,  or  65  meters  at  mid¬ 
latitude. 

Above  about  80  to  90  kilometers,  one  must  be  aware  of  the 
distinction  between  kinetic  and  molecular-scale  temperature. 

Earlier  models  used  the  latter,  while  later  models  generally  use 
the  former.  The  error  increases  monotonically  with  increased 
altitude  and  can  be  as  high  as  1500°K  in  the  thermosphere. 

As  mentioned  previously,  many  reference  atmospheres  are 
moist  models,  and  care  must  be  used  in  the  distinction  between 
kinetic  and  virtual  temperature  at  lower  levels  (below  about  10  km) . 
Although  the  difference  can  be  as  high  as  6°K,  for  practical 
purposes,  it  seldom  exceeds  3°K.  It  is  important  to  note,  for 
example,  that  the  tabular  values  given  in  the  U.S.  Standard 
Atmosphere  Supplements,  1966,  are  virtual  temperatures  which  are 
derived  from  relative  humidity  values  given  in  Table  E-l. 
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The  final  major  source  of  error  is  in  the  algorithms  used 
to  reproduce  standard  atmospheres  for  computer  applications.  Of 
course,  care  must  be  taken  to  account  for  the  different  heights 
and  temperatures  discussed  above.  There  are  two  basic  types  of 
algorithms  most  often  used.  The  first,  application  of  che 
integral  form  of  the  hydrostatic  equation  to  the  temperature 
profile,  is  generally  the  most  accurate,  but  is  also  the  most 
time  consuming  unless  data  is  requested  in  order  of  increasing 
altitude.  The  second  general  method  involves  a  table  lookup  or 
interpolation  on  two  parameters  (e.g.  temperature  and  pressure) 
and  calculation  of  the  third  from  the  perfect  gas  law.  To 
accurately  interpolate  on  pressure  or  density,  tabular  values 
must  be  available  at  altitude  intervals  of  from  less  than  about 
250  meters  near  the  surface  to  20  kilometers  in  the  thermosphere. 
The  obvious  disadvantage  of  this  technique  is  the  large  data  arrays 
required.  It  is,  however,  generally  faster  than  the  first  method. 
Round  off  error  is  generally  not  a  problem  on  most  machines  unless 
some  kind  of  polynomial  curve  fit  is  employed,  or  the  first 
method  is  employed  to  great  heights  at  relatively  small  intervals. 
In  such  cases,  double  precision  variables  should  be  used, 
d .  Meteorological  Constants 

Table  E-2  is  a  summary  of  the  primary  constants  adopted  for 
the  U.S.  Standard  Atmosphere,  1962.  They  are  presented  here  as 
standards  for  meteorological  calculations,  to  be  used  when  reducing 
data  from  the  standard  to  conditions  compatible  with  other  data. 
Some  of  them  apply,  strictly,  only  to  4  5°N  latitude.  Conditions 
at  other  latitudes  can  be  obtained  from  supplementary  atmospheres 
such  as  the  U.S.  Standard  Atmosphere  Supplements,  1966. 

Supplementary  or  derived  constants  are  given  in  Table  E-3. 
Abundances  of  the  most  common  atmospheric  constituents  of  dry  air 
are  given  in  Table  E-4.  Other,  less  used,  constants  and  equations 
are  given  in  the  U.S.  Standard  Atmosphere,  1962  and  the  U.S. 
Standard  Atmosphere  Supplements,  1966, 
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Figure  E-l 

Temperature-height  protiles  of  the  U.S.  Standard  Atmosphere,  1962 
compared  with  the  ARDC  1956  and  1959  and  the  ICAO  1952  and  1954.  The  COESA 
1958  agrees  closely  with  the  ARDC  1956.  There  are  additional  small  differ¬ 
ences  due  to  changes  in  some  physical  constants.  Pressure  and  density  scales 
refer  to  the  1962  standard.  (After  Valley,  1956) 


TABLE  e-1.  MOISTURE  PROPERTIES  OF  THE  1966 
U.S.  STANDARD  ATMOSPHERE  SUPPLEMENTS  (after 
Valley,  1956) 


Altitude 

(km) 

Temperature 

Virtual 

Relative 

Geom 

Geop 

( °K) 

Temp  (°K) 

Humidity,  % 

Tropical 

(15°  N) 

0.000 

0.000 

299.65 

302.588 

75 

1.002 

1.000 

293.65 

295.893 

75 

2.005 

2.000 

287.65 

289.336 

75 

2.256 

2.250 

286.15 

287.717 

75 

2.507 

2.500 

286.95 

287.743 

35 

4.012 

4.000 

276.90 

277.363 

35 

6.020 

6.000 

263.50 

263.709 

35 

8. 029 

8.000 

250.10 

250.172 

30 

10.039 

10.000 

236.70 

236.717 

20 

Temp  (°K) 

Virtual  Rel 

Jan  July 

Temp 

(  K)  Hum  (%) 

Jan 

July  Dan  July 

Subtropic  (30°N) 


0.000 

0.000 

287.15 

301.15 

288.519 

304.583 

80 

80 

1.002 

1.000 

284.15 

293.65 

285.244 

295.580 

70 

65 

2.003 

2.000 

281.15 

288.15 

281.862 

289.536 

50 

60 

3.006 

3.000 

274.65 

282.65 

275.098 

283.716 

45 

60 

4.008 

4.000 

268.15 

277.15 

268.389 

277.823 

35 

50 

6.014 

6.000 

255.15 

266.15 

255.239 

266.445 

30 

40 

8.021 

8.000 

242.15 

252.15 

242.185 

252.266 

30 

40 

10.030 

10.000 

299.15 

238.15 

229.162 

238.179 

30 

30 

Midlatitude 

(45°N) 

0.000 

0.000 

272.15 

294.15 

272.594 

296.216 

77 

75 

1.000 

1.000 

268.65 

289.65 

268.998 

291.142 

70 

65 

2.001 

2.000 

265.15 

285.15 

265.427 

286.192 

65 

55 

3.001 

3.000 

261.65 

279.15 

261.850 

279.777 

55 

45 

4.003 

4.000 

255.65 

273.15 

255.774 

273.552 

50 

40 

6.005 

6.000 

24  3 . 65 

261.15 

243.698 

261.299 

45 

30 

8.010 

8.000 

233 .65 

248.15 

231.664 

248.211 

35 

30 

10.016 

10.000 

219.65 

235.15 

219.654 

235.172 

30 

30 

Subarctic 

(60°N) 

0.000 

0.000 

257.15 

287.15 

257.285 

288.449 

80 

7  5 

0.999 

1.000 

259.15 

281.75 

259.311 

282.685 

70 

70 

1.998 

2.000 

255.95 

276.35 

256.089 

277.062 

70 

70 

2.998 

3.000 

252.75 

270.95 

252.861 

271.447 

65 

65 

3.497 

3.500 

251.15 

268.25 

251.245 

# 

60 

— 

3.997 

4.000 

247.75 

265.55 

247.824 

265.889 

60 

60 

4.998 

5.000 

240.95 

260.15 

# 

260.376 

— 

55 

5.998 

6.000 

234.15 

253.15 

234.170 

253.277 

50 

50 

8.000 

8.000 

220.55 

239.15 

220.550 

239.185 

40 

40 

10.003 

10.000 

225.15 

225.155 

30 

h 


t 


TABLE  E- 1 .  ( CONCLUDED) 


Altitude 

Geom 


(km) 

Geop 


Temp  (  K) 
Jan  July 


Virtual 
Temp  (°K) 
Jan  Jul5 


Arctic  (75°N) 


Rel 

Hum  ( % ) 
Jan  July 


0.000 

0.000 

249.15 

278.15 

249.216 

278.924 

80 

85 

0.998 

1.000 

252.15 

275.55 

252.231 

276.187 

65 

75 

1.497 

1.500 

253.65 

274.25 

253.741 

# 

60 

— 

1.996 

2.000 

250.90 

272.95 

250.976 

273.463 

60 

65 

2.495 

2.500 

248.15 

271.65 

# 

272.144 

— 

65 

2.995 

3.000 

245.40 

268.40 

245.448 

# 

55 

— 

3.994 

4.000 

239.90 

261.90 

239.929 

262.131 

50 

55 

5.992 

6.000 

228.90 

248.90 

228.911 

248.978 

45 

45 

7.992 

8.000 

217.90 

235.90 

217.90 

235.922 

40 

35 

9.493 

9.500 

266.15 

226.158 

— 

30 

9.993 

10.000 

226.65 

226.656 

— 

20 

#  Not  a  virtual  temperature  breakpoint, 
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TABLE  E-2 

PRIMARY  METEOROLOGICAL  CONSTANTS  USED  IN 
DEVELOPING  1962  U.S.  STANDARD  ATMOSPHERE 


Quantity 

Symbol 

Value 

Error  ^ 

Units 

Sea-level  pressure 

(45°N) 

po 

1.01325  x 

105 

defined 

N-m'2  (2) 

Sea-level  density 

(45°N) 

1.22500  x 

10"7 

5 

,  -3 

kg*m 

Sea-level  temperature  T^ 
(45°N) 

288.15  (15 

.00) 

defined 

K  ( °C) 

Sea-level  accelera¬ 
tion  ( 45°N) 

^0 

9.80665 

defined 

-2 

m*  s 

Ice-point  temperature  T. 

273.15 

defined 

°K 

Triple-point  tempera 
ture  (H20) 

—  ““ 

273.16 

8 

defined 

°K 

Mean  collision 

diameter  (air) 

3.65  x  10" 

1 

m 

Avogadro's  Number 

na 

6.02257  x 

io26 

defined 

(3) 

kmole  1 

Universal  Gas 

Constant 

R 

8.31432  x 

io3 4 

4 

J-kmole”1^"1 

Molecular  weight  of 

water 

M 

w 

18.0153 

1 

i  ,  i  -1  (4) 
kg -kmole 

l 


t 

l 

h  ! 


(1)  Value  is  plus  or  minus  the  last  digit  given 

(2)  Ordinarily  given  as  1013.25  mb  (millibars) 

(3)  Defined  for  purposes  of  these  calculations,  differs  from  SI 
value  by  .0004 

(4)  Based  on  C"^2  =  12.0000  kg*kmole  ^ 


i 

I:, 
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TABLE  E-3 


DERIVED  METEOROLOGICAL  CONSTANTS  USED  IN 
DEVELOPING  1962  U.S.  STANDARD  ATMOSPHERE 


Quantity 

Symbol 

Value 

Error 

Units 

Molecular  weight 

dry  air 

of 

Md 

28.9644 (2) 

2 

defined 

kg*kmole  ^ 

-1„  -i 

Gas  constant  for 

air 

dry 

Rd 

2.87053  >: 

10 

defined 

J-kg  •  K 

Specific  heat  of 

air,  constant 

pressure 

dry 

CP 

10.04686  >: 

:  102 

(1) 

J-kg_1-0K-1 

Specific  heat  of 

air,  constant 

volume 

dry 

C 

V 

7.17633  >: 

102 

(1) 

J-kg'^K"1 

Ratio  of  specific 

heats  (C  /C  ) 

p  V 

If 

1.4 

defined 

(dimension¬ 

less) 

(1)  Value  is  plus  or  minus  the  last  digit  given 

(2)  Based  on  C^  =  12.0000  kg*kmole  ^ 
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TABLE  E-4 


NORMAL  COMPOSITION  OF  CLEAN,  PRY  AIR  NEAR  SEA  LEVEL 


Gas/ 

Formula 

By  Volume 

By  Weight 

Molecular 
Weight (1) 

Nitrogen  (N2) 

78.08410.004% 

75.520210.004% 

28.0134 

Oxygen  ( O^ ) 

20.910*0.002% 

23.140410.002% 

31.9988 

Argon  (A) 

0.93410.001% 

1.28810.001% 

39.948 

Carbon  Dioxide 

(C02) 

*0.03310.001% 

*0.05010.002% 

44.00995 

Neon  (Ne) 

18 . 1810 . 04ppm 

12 . 6710 . 04ppm 

20.183 

Helium  (He) 

5.2410. 04ppm 

0 .72410 . 006ppm 

4.0026 

Methane  (CH^) 

*2  ppm 

*lppm 

16.04303 

Hydrogen  (H2) 

0.510. 05ppm 

0 . 0310 . 003ppm 

2.01594 

Nitrous  Oxide 

«20) 

*  0 . 51 0 . lppm 

*0.810.2  ppm 

44.0128 

Xenon  (Xe) 

0 . 08710 . OOlppm  0 . 394 10 . 005ppm 

131.30 

Ozone  (03) 

*0  to  0.07ppm 

*0  to  O.lppm 

47.9982 

Sulfur  Dioxide 

(S02) 

*0  to  10. ppm 

*0  to  20. ppm 

64.0628 

Nitrogen  Dioxide  (N02) 

*0  to  0.02ppm 

*0  to  0.03ppm 

46.0055 

Iodine  (I_) 

*0  to  O.Olppm 

*0  to  0.09ppm 

253.8088 

(1)  based  on  C ^  =  12.0000  kg*kmole  1 

(*)  somewhat  variable  due  to  pollution 
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ASD  TR-79-5056 


LOCATION 

ASD/WE 

Wright -Patterson  AFB,  OH  45433 
AFWAL/AATS 

Wright-Patterson  AFB,  OH  45433 
ASD/PAM 

Wright-Patterson  AFB,  OH  45433 
AFEWC/EST  ' 

San  Antonio,  TX  78243 
AFSC/IN 

Andrews  AFB,  Wash  D.C.  20334 

Air  Univ  Library 
Maxwell  AFB,  AL  36112 

J  DTIC 

*  Cameron  Station,  Alexandria,  VA 

22314 

USAF/SAMID 

Washington,  D.C.  20330 
ASD/ENFTC 

Wright-Patterson  AFB,  OH  45433 

USAFETAC/TS 
Scott  AFB,  IL  62225 

USAFETAC/DN 
Scott  AFB,  IL  62225 

AFGWC/TS 

Offutt  AFB,  NE  68113 
2WS/DR 

Andrews  AFB,  Wash  D.C.  20334 
FTD/WE 

Wright-Patterson  AFB,  OH  45433 
ASD/ADH 

Wright-Patterson  AFB,  OH  45433 
AFGL/XOP 

Hanscom  AFB,  MA  01730 


NO.  COPIES 
18 

2 

1 

3 

1 

1 

12 


1 

25 

1 

1 

1 

1 

1 

2 

1 


AFWAL/AARI-3 

Wright-Patterson  AFB,  OH  45433 
AFWAL/KEA 

Wright-Patterson  AFB,  OH  45433 
AFWAL/WEF 

Wright-Patterson  AFB,  OH  45433 

Det  1,  2WS  (1RF) 
Wright-Patterson  AFB,  OH  45433 

ASD/XRH 

Wright-Patterson  AFB,  OH  45433 
ASD/YWE 

Wright-Patterson  AFB,  OH  45433 
ASD/SD28D 

Wright-Patterson  AFB,  OH  45433 
ASD/SD30E 

Wright-Patterson  AFB,  OH  45433 

ASD/SD30EF  (Mr.  Simerlink) 
Wright-Patterson  AFB,  OH  45433 

O.E.  Smith 
CODE  ES-81 
NASA 

Marshall  Space  Flight  Center 
Huntsville,  AL  35812 

AFWAL/FIEE 
ATTN:  Dr.  Mayer 
Wright-Patterson  AFB,  OH  45433 

AFWAL/FIEE-A 
ATTN:  Mr.  Warner 
Wright-Patterson  AFB,  OH  45433 

ASD/ENO 

Wright-Patterson  AFB,  OH  45433 
ASD/YF/EC 

Wright-Patterson  AFB,  OH  45433 


